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Abstract 

Weak scale super symmetry is often said to be fine-tuned, especially if tfie matter content 
is minimal. This is not true if there is a large A term for the top squarks. We present a 
systematic study on fine-tuning in minimal super symmetric theories and identify low energy 
spectra that do not lead to severe fine-tuning. Characteristic features of these spectra are: 
a large A term for the top squarks, small top squark masses, moderately large tan/3, and 
a small /i parameter. There are classes of theories leading to these features, which are 
discussed. In one class, which allows a complete elimination of fine-tuning, the Higgsinos 
are the lightest among all the superpartners of the standard model particles, leading to 
three nearly degenerate neutralino/chargino states. This gives interesting signals at the 
LHC — the dilepton invariant mass distribution has a very small endpoint and shows a 
particular shape determined by the Higgsino nature of the two lightest neutralinos. We 
demonstrate that these signals are indeed useful in realistic analyses by performing Monte 
Carlo simulations, including detector simulations and background estimations. We also 
present a method that allows the determination of all the relevant superparticle masses 
without using input from particular models, despite the limited kinematical information 
due to short cascades. This allows us to test various possible models, which is demonstrated 
in the case of a model with mixed moduli-anomaly mediation. We also give a simple 
derivation of special renormalization group properties associated with moduli mediated 
supersymmetry breaking, which are relevant in a model without fine-tuning. 



1 Introduction 



What is the physics at the TeV scale and how can we test it? These questions become more 
pressing as we approach the LHC era, which will start within two years. It is extremely important 
now to consider what we expect to see at these energies, especially because the LHC is a hadron 
collider experiment, in which relations between experimental data and the underlying theory are 
not so simple. Knowing what we are looking for would certainly help to identify the physics at 
the TeV scale and may even be necessary, as the determination of the TeV physics at the LHC 
and other experiments will most likely take the form of a slow elimination process. 

There are already several hints on possible physics at the TeV scale. They come from combin- 
ing a theoretical criterion of naturalness and precision measurements of electroweak observables 
and rare flavor- and CP- violating processes. Among these, the combination of naturalness and 
the electroweak data seems to give the most unambiguous hint, because these constraints cannot 
be satisfied simply by imposing approximate symmetries already present in the standard-model 
gauge and matter sector. Interpreted naively, the precision electroweak data suggest the existence 
of a light Higgs boson, with the contributions to the electroweak observables from other physics 
highly suppressed [1]. Naturalness then implies that there must be a new weakly-interacting 
physics at a TeV scale or below, which cuts off quadratically divergent contributions to the 
Higgs mass-squared parameter arising from loops of the standard model particles. 

Weak scale supcrsymmetry is an ideal candidate for the new physics. Loops of supcrparticlcs 
cancel the quadratic divergences from those of the standard model particles. The new inter- 
actions for the superparticles are necessarily weak, as they are related to the standard model 
interactions by supcrsymmetry. Moreover, the mass of the lightest Higgs boson is predicted to 
be small, Mniggs ^ 200 GeV in most (even extended) theories, which is very much consistent 
with the precision electroweak data. With the introduction of R parity and the assumption of 
flavor universality for the superparticle masses, weak scale supersymmetry can provide a fully 
consistent framework for physics of electroweak symmetry breaking. 

Postulating weak scale supersymmetry alone, however, does not much narrow down signa- 
tures at the LHC. Depending on the relative sizes for the soft supersymmetry breaking pa- 
rameters, one can have drastically different signatures at the LHC. The number of relatively 
model-independent signals is also small, making it difficult to discriminate supersymmetry from 
other TeV-scale physics. A generic signal of weak scale supersymmetry is large missing trans- 
verse energy in association with jets and/or isolated leptons. Such a signal, however, arises in 
almost any theory where the lightest TeV-scale particle is stable and neutral, which is suggested 
by the existence of the dark matter of the universe. It is then an important task to narrow down 
the parameter space of weak scale supersymmetry further and to perform a detailed study of 
the LHC signals there. One of the important goals of such a study is to identify generic signals 
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associated with a particular parameter region so that the non-observation of those signals will 
allow us to exclude the region. 

How should we choose regions among the huge parameter space of weak scale supersymmetry? 
An obvious way is to assume a particular supersymmetry breaking model, such as the minimal 
supergravity or gauge mediation model, based on the simplicity of the model. This selects a 
slice in the parameter space of soft supersymmetry breaking parameters, which depends only 
on a few free parameters. These studies have been performed by many authors, for example, 
in [2 - 8]. In this paper we take a different criterion to choose the region. We use the hint 
from naturalness to the maximal amount and consider what are generic implications of it on the 
spectrum of superparticles and on LHC signals. Fortunately, or unfortunately, generic parameter 
regions of weak scale supersymmetry satisfying existing experimental constraints do not lead 
to the correct scale for electroweak symmetry breaking without significant fine-tuning. This 
information, therefore, can be used to constrain the parameter space of soft supersymmetry 
breaking parameters and thus to narrow down possible signatures at the LHC. Of course, the 
input from naturalness alone does not lead to unambiguous signatures in a wide variety of 
possible supersymmetric theories. In this paper we focus our attention to the case where the 
matter content at the weak scale is minimal, i.e. given by that of the minimal supersymmetric 
standard model (MSSM). 

What are generic implications of fine-tuning on the spectrum of superparticles in a theory 
with the minimal matter content? As discussed in Ref. [9], the fine-tuning problem of minimal 
supersymmetry can be solved without extending its matter content if the trilinear scalar term 
{A term) for the top squarks is large and the holomorphic supersymmetry-breaking term {fiB 
term) for the Higgs doublets is small. This allows us to evade the LEP II bound on the Higgs 
boson mass with relatively small superparticle, specifically top squark, masses. Then, if soft 
supersymmetry breaking parameters are generated (effectively) at low energies, the sensitivity 
of the electroweak scale to the fundamental parameters of the theory can be very mild. One 
of the consequences of such a scenario is that the top squarks arc relatively light and have a 
large mass splitting between the light and heavy ones. Another important consequence is that 
the Higgsinos are rather hght, with the masses smaller than about 190 GeV (270 GeV) for fine- 
tuning better than fa 20% (10%), which is because naturalness requires any contribution to the 
Higgs boson squared mass to be small, including the supersymmetric contribution. We argue 
that these features are robust and appear quite generically in a minimal supersymmetric theory 
without significant fine-tuning. 

This argument provides a strong motivation to consider the case in which the lightest neutral 
Higgsino is the lightest supersymmetric particle (LSP). While it is not a necessary consequence of 
solving the supersymmetric fine-tuning problem, the Higgsino LSP in fact arises in a large class 
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of theories in which the pattern of soft supersymmetry breaking parameters described above is 
naturally obtained. It is, therefore, quite important to perform an LHC study for the case of the 
Higgsino LSP and identify possible signatures. An important consequence of the Higgsino LSP 
scenario is that there are three nearly degenerate neutralino/chargino states, Xi, X2 xt^ with 
Xi being the LSP. We show that this structure can give interesting signatures at the LHC in the 
dilcpton invariant mass distribution arising from the decay X2 ~^ Xi^~^^~- We discuss in what 
sense these arc characteristic signatures of the Higgsino LSP, and under what circumstances the 
signals can be used in realistic analyses. 

To demonstrate the usefulness of the signatures, we need to choose specific parameter points 
and perform Monte Carlo simulations, including detector simulations and standard model back- 
ground. We do this in the model discussed in Refs. [9 - 11], where the desired pattern of the 
soft supersymmetry breaking masses, a large A term and a small /iB term, are obtained while 
evading the existing experimental constraints such as the one from b ^ s^. We show that the 
dilepton invariant mass distribution from X2 ~^ Xi is indeed useful to test the Higgsino 
nature of the LSP and to extract the information on a small mass difference between xi and 
X2- We also show that important parameters of the model, the overall mass scale and the // 
parameter, are determined by various other endpoint analyses. In fact, we show that these pa- 
rameters are overconstrained, so that we can test some of the model predictions. We perform 
these analyses for an integrated luminosity of 30 fb~^, but essentially the same conclusion is ob- 
tained with 10 fb~^. The technique presented here can also be used in a larger class of theories 
having similar superparticle spectra. 

The organization of the paper is as follows. In section 2, we present a systematic study on 
naturalness in general supersymmetric theories, especially focusing on the case where the matter 
content at the weak scale is minimal. We give general criteria that natural supersymmetric 
models with the minimal matter content must satisfy, and present characteristic patterns for the 
superparticle spectrum arising from these models. In section 3, we discuss LHC signals of the 
Higgsino LSP scenario, which naturally arises in a class of models that do not suffer from fine- 
tuning. We find that a combination of the endpoint and the shape of the dilepton invariant mass 
distribution provides a powerful tool to test the scenario. In section 4, we perform Monte Carlo 
simulations to demonstrate that these signals are indeed useful in a realistic situation. We also 
present an analysis that allows us to determine the masses of the gluino, squarks, and the two 
lightest neutralinos in the Higgsino LSP scenario, without relying on details of the underlying 
model. We illustrate that these information can be used to test and/or discriminate between 
possible models. Discussion and conclusions are given in section 5. In the Appendix, we give a 
simple derivation of special renormalization group properties of moduli mediated supersymmetry 
breaking models, which are relevant in the model studied in section 4. 
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2 Supersymmetry and Naturalness 



One of the principal motivations for weak scale supersymmetry is to provide a solution to the 
naturalness problem of the standard model. This has, however, been put in a subtle position after 
non-discovery of both superparticles and a light Higgs boson at LEP II. In a generic parameter 
region motivated by simple supersymmetry breaking models, fine-tuning of order a few percent 
is required to reproduce the correct scale for electroweak symmetry breaking while evading the 
constraints from LEP II. This problem, called the supersymmetric fine-tuning problem, has 
attracted much attention recently, and several solutions have been proposed, e.g., in [9, 10, 12 - 
22]. In this section we reconsider the problem and see what are generic implications of it, 
especially in the context of theories with the minimal matter content. One of our emphases 
here is on the fact that the supersymmetric fine-tuning problem may simply be a problem of the 
supersymmetry breaking mechanism and not necessarily that of minimal supersymmetry itself. 

2.1 Large At and small y^B in minimal supersymmetry 

We begin our discussion by considering fine-tuning in the Higgs potential in general weak scale 
supersymmetric theories. Let h be the Higgs field whose vacuum expectation value (VEV) 
breaks the electroweak symmetry. In minimal supersymmetry, /i is a linear combination of the 
two Higgs doublets, B.u and B.^,- The potential for h is given by 

V = ml\h\^^^-l\h\\ (1) 

where m| is negative and \h is positive. By minimizing it, wc obtain = {h^ — —2m\/\h and 
-^Higgs — ^hV^, where Mniggs is the mass of the physical Higgs boson, so that 

= -ml (2) 

We thus find that |m^| cannot be large for a light Higgs boson: |m^|^/^ ^ 140 GeV (90 GeV) 
for Mffiggs ^ 200 GeV (130 GeV). 

What is ml in supersymmetric theories? For moderately large tan/3 = {Hy)/{H(i), e.g. 
tan /? ^ 2, ml can be written as 

ml = l/^l^ + mifjtree + "^kl^ad, (3) 

where /i is the supersymmetric mass for the Higgs doublets, and mjjjtree and mj^J^ad represent 
the tree-level and radiative contributions to the soft supersymmetry-breaking mass squared for 
H^. The dominant contribution to m|^^|rad arises from top-stop loop: 

2 I 3|/f / 2 I 2 I I /I |2\ 1 ( -^mcss \ /^\ 

%Jrad--^(mQ3+m^3 + |A| )ln — — 1, (4) 



m. 
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where yt is the top Yukawa couphng, rriq^ and m^^ soft supersymmetry breaking masses for the 
third-generation doublet quark, Q3, and singlet up-type quark, f/3, and At the trilinear scalar 
interaction parameter for the top squarks (our definition for the A parameters is such that a 
scalar trilinear coupling is given by the product of the Yukawa coupling and the A parameter, 
e.g., C = —ytAtq^u^Hu + h.c). The quantity Mmcss represents the scale at which squark and 
slepton masses are generated, and mi the scale of the top squark masses determined by toq^, 
mfj^ and Af. Note that Mmcss can be an effective scale different from the true scale of scalar 
mass generation in a case that the theory possesses special relations among various parameters. 

For fine-tuning to be absent, each term in the right-hand-side of Eq. (3) should not be much 
larger than the left-hand-side, which is related to the physical Higgs boson mass by Eq. (2). Let 
us first consider m^^|rad- The amount of fine-tuning from this term is given by ^Higgs/2^//„|rad, 
so that requiring the absence of fine-tuning worse than leads to the condition 

where we have set TOq^ ~ mfj,^ ~ m? for simplicity, and x = \At\/mi. This has the following 
implication on the properties of the supersymmetry breaking sector [18]. Unless Mmcss is ex- 
tremely small, e.g. Mmcss ^ 10 TeV, the absence of fine-tuning, defined by A^^ > 20%, requires 
mi ^ 700 GeV, where we have used Mniggs ^ 200 GeV as suggested by the precision electroweak 
data. This implies that a naive low-scale mediation model, leading to the "minimal gauge me- 
diated mass relation" m?/m? ^ ot/gf, is unlikely to solve the fine-tuning problem because it 
gives too light right-handed sleptons. For a lighter Higgs boson, we obtain severer bounds on 
mf for Mniggs — 140 GeV, for example, we find mi ^ 700 GeV even for ~ 10%. Note that 
the condition of Eq. (5) applies independently of any other considerations.''^ 

Light top squarks suggested by Eq. (5) leads to a tension with the LEP II bound on the Higgs 
boson mass, Mniggs ^ 114.4 GeV [23], since in the MSSM having Mniggs larger than the Z boson 
mass, mz, requires radiative corrections arising from top-stop loop, which grow with the top 
squark masses [24] . A simple way to avoid the conflict is to introduce an additional contribution 
to the Higgs boson mass other than that in the MSSM. An example of such theories can be found 
in Ref . [18] , where the required properties for the supersymmetry breaking sector are realized by 
strong gauge dynamics breaking supersymmetry at a scale of (10 ~ 100) TeV. What if we do not 
introduce any other contribution to Mniggs than that in the MSSM? In this case it is unlikely 

^In the case that hi(Mnicss/?72() is large, for example in gravity mediated models, the expression in Eq. (4) is 
not reliable and we should sum up the leading logarithms using renormalization group equations. This case will 
be addressed in the next subsection. 
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Figure 1: Minimal values of = (TOQ^y^'^ = {mfj^Y^'^ giving Mniggs ^ 114.4 GeV as a function 
of At/ mi. The other parameters are fixed to be 500 GeV for the gaugino and sfermion masses 
other than m^^ and Tnfj^, (500 GeV)(Af/mj) for the A parameters other than At, tan/? = 15, 
H = +170 GeV, and = 250 GeV. 

that Mniggs can be larger than 130 GeV, so Eq. (5) leads to a severer bound 

mf S (450 GeV)^^ (^^^ ^ (6) 

While this bound is strong, it still leaves a room for evading the LEP II constraint on Mniggs 
without introducing severe fine-tuning. As discussed in [9], this happens if At is large, tan/3 is 
(moderately) large, and Mmess is small. In particular, it is crucial to have large At, compared 
with rrit, to evade the Higgs boson mass bound while keeping modest. 

In Fig. 1, we plot minimal values of m^- = ("^Qg)^^^ = (nT'i/.J^^'^, ""^tlmin, that give M^ggs ^ 
114.4 GeV as a function of At/ mi. The other parameters are fixed to be 500 GeV for the gaugino 
and sfermion masses other than m^^ and m\j^, (500 GeV){At/mi) for the A parameters other 
than At, tan/3 = 15, = +170 GeV, and mA = 250 GeV, where is the mass of the pseudo- 
scalar Higgs boson. The dependence of the results on the fixed parameters is not significant, 
as long as tan/5 is sufficiently large, e.g., tan/5 ^ 10. In the figure we plot mf|min for three 
different values of the top quark mass, mt = 166.9, 172.7 and 178.5 GeV, corresponding to the 
central value and the 2a range of the the latest experimental data: mt = 172.7 ± 2.9 GeV [25]. 
The calculation here has been performed using FeynHiggs 2.2 [26] (for earlier analyses for the 
Higgs boson mass in the MSSM, see e.g. [27]). If we instead use a code based on the pure DR 
scheme, such as SuSpect 2.3 [28], we obtain slightly different values of 
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the differences are of order 50 GeV but for \At/mf\ -C 1 the DR scheme calculation can give 
m^lmin larger than that in the figure by about 200 GeV. These difi^erences give an estimate for 
the size of higher order corrections. Throughout the paper, our sign convention for n and the 
soft supersymmetry breaking parameters follows that of SUSY Lcs Houches Accord [29]. 

The figure clearly shows that in order to have light top squarks suggested by Eq. (6) the 
existence of a substantial At term (|At/mf| ^ 1) is required. For nit = 172.7 GeV, the existence 
of At with lAt/m^l ^ 1 can allow as small as (200 ~ 400) GeV while should be larger 
than ^ 1.2 TeV for \At/mt\ <^ 1. Another important point is that for |At/m^| ~ 2, which 
gives the minimal value of mtlnim, the sensitivity of mj|min to the value of rrit is mild, while for 
\At/mt\ <^ 1 the sensitivity is huge. This implies, for example, that if turns out to be smaller 
than ~ 170 GeV theories with \At/mt\ -C 1 at the weak scale will pretty much be "excluded." 

We conclude that to have natural electroweak symmetry breaking in super symmetric theories 
with the minimal, i.e. MSSM, matter content, the existence of a substantial At term at the weak 
scale is crucial. Another important ingredient is a moderately large tan (3, e.g. tan /3 ^ 5, to have 
a sufficiently large tree-level Higgs boson mass, which requires the holomorphic supersymmetry 
breaking mass squared for the Higgs doublets, the /iB term, to be (significantly) smaller than 
2|/ip + rnjj^ + Tn?Hd- ^St^ci, these ingredients dominantly control the amount of fine-tuning in 
almost any theory with the MSSM matter content. To demonstrate this, we will now analyze 
the situations in the case of high scale supersymmetry breaking and in gauge mediation models 
from the viewpoint of the size of At at the weak scale. For earlier analyses of fine-tuning in these 
models, see e.g. [30]. 

2.2 High scale supersymmetry breaking and gauge mediation 

Let us first consider the minimal supergravity (mSUGRA) scenario [31]. We start by considering 
the constrained mSUGRA, in which the soft supersymmetry breaking parameters are specified 
by the universal gaugino mass M1/2, universal scalar mass squared ttIq, universal A term Aq, 
and the iiB term at the unification scale, Munif ~ 10^^ GeV. While this scenario is sometimes 
criticized due to a lack of a strong theoretical motivation, it is not so bad in term of fine-tuning, 
compared with other models such as gauge mediation models. This is because we can obtain a 
reasonable size of At /mi at the weak scale, so that the top squark masses can be made smaller 
compared with the models giving smaller values of At/ mi at the weak scale. 

In Fig. 2 we plot values of the fine-tuning parameter A^^ in the constrained mSUGRA for 
two different choices of Aq. = and Aq = — 3|mo|. The parameter is defined by the 
fractional sensitivities of the electroweak VEV, f ^ 174 GeV, to the fundamental parameters 
of the theory, with generic sensitivities of v to the parameters appropriately corrected [32]. We 
plot the contours of A^^ on the mQ-Mi/2 plane for /x > 0, where mo = sgn(mQ)|mQ|^/^. The 
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Figure 2: Contours of A^^ on the mo-Mi/2 plane for the constrained mSUGRA with = 
(left) and Aq = — 3|mo| (right). The sign of /z is chosen to be positive. The constraints from 
direct superparticle search, the Higgs boson mass bound, and the stau LSP are also shown. 

values of /i and at Munif are determined by v and tan/3, and we take tan/? = 15. We 
find that for Aq = the fine-tuning is worse than 2%, while for Aq = — 3|mo| it can be as 
mild as 5% for M1/2 — 150 GeV and ttiq ~ (200 GeV)^. This can be understood as follows. 
For y4o = 0, renormalization group equations give low-energy values for At and m^- that satisfy 
At /nil ~ — 1- While this value of |At/mj| is not totally negligible, it is still not large enough to 
give Mniggs ^ 114.4 GeV with top squark masses smaller than about 600 GeV (see Fig. 1). This 
gives a high sensitivity of v to yt (the top-stop contribution to rrijj^), leading to ^ 2%. The 
situation can be made better by introducing non- vanishing Aq at Munif. While the sensitivity of 
low-energy At to Aq is rather weak, Aq = — 3|mo| can give a low-energy value of At/mi about 
—1.8, which allows mi as small as ~ 250 GeV to evade the Higgs boson mass bound, and thus 
as large as 5%. Here mi is defined by mi = (mg^m^^)^/^. In fact, larger values of Aq do 
not help in reducing fine-tuning because of a shrinking of the phenomenologically acceptable 
parameter region, and we obtain A~^|inax ~ 5% in the constrained mSUGRA. 

In the case of the constrained mSUGRA described above, A~^ is determined by the sensitivity 
of V to Ut and /i, which implies that the dominant source of fine-tuning comes from the sensitivity 
of m^j^ to the top-stop loop contribution. We can make this sensitivity weaker by deviating from 
the constrained mSUGRA. A simple way of doing this is to make m^j^ and differ from m^ 
at Munif. Practically, this implies that we can take low-energy values of ^ and as free 
parameters. Then, for certain values of and m^i, which corresponds to choosing certain values 
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Figure 3: Contours of A^^ on the mo-Mi/2 plane in the mSUGRA model with = 190 GeV 
and mA = 250 GeV fixed at the weak scale and Aq = — 3|mo| at Munif (left). Contours of 
on the mg-Mi/2 plane with > and m? = (500 GeV)^, Aq = —\mg\ and 171^ = m^^ = m^^ = 
(100 GeV)^ at Munif, where m? and m? are the squark and slepton masses, respectively (right). 
In both cases tan/5 = 15. 

of rriff^ and m|^^ at Munif, we find that the sensitivity of v to yt can be made weaker due 
to renormalization group properties of the soft supersymmetry breaking parameters. This is 
illustrated in the left panel of Fig. 3, where we plot the contours of on the mo-Mi/2 plane, 
with fi = 190 GeV, = 250 GeV and Aq = — 3|mo|. We find that fine-tuning can be made as 
mild as A^^ ^ 8%. A similar reduction of tuning can also occur in the region with m? < at 
Munif, which is allowed if we violate the universality of M1/2 and/or ttLq to avoid the slepton mass 
bound. This is illustrated in the right panel of Fig. 3 for yU > 0, where we plot A~^ as a function 
of the gaugino mass, M1/2, and the squark mass, rrig = sgn(m?)|m?|^/^, at the unification scale 
Munif. The other parameters are chosen as m? = (500 GeV)^ for the sleptons, Aq = — \mq\, and 
mjj^ = mjj^ = (100 GeV)^ at Munif. We find that fine-tuning can be as mild as A"^ ^ 8%.^ 

How much can we reduce fine-tuning in a theory with high scale supersymmetry breaking? It 
is possible, after all, that the soft supersymmetry breaking parameters are not calculable (easily) 
if they are generated through physics at the gravitational scale. Suppose, for example, that grand 
unification is realized in five dimensions and supersymmetry is broken on a brane on which the 

^The region with to? < has been discussed recently in [22] in the context of finding a relation among soft 
supersymmetry breaking parameters which reduces fine-tuning. Our approach here is different: we do not assume 
any special relations among the supersymmetry breaking masses, e.g., between the gaugino and squark masses. 
We then do not find a region with better than « 10%. 
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active gauge group is only SU (3)c x SU{2)l x U{1)y [33]. Suppose also that all the gauge, matter 
and Higgs fields propagate in the bulk, so that they all feel supersymmetry breaking through the 
operators [ZW^Wale^ , [{Z -\- Z'^)^^^]g4, and [Z"'"Z$''"$]g4, where Z is the supersymmetry breaking 
field, Wa the gauge field-strength superfields, and $ the matter and Higgs chiral superfields. 
(The n and fiB terms can also be generated through [Z'^ HuH^-\-h.c.]g4 and [Z^^ ZHuHd + h.c.]g4.) 
Then, if there is a flavor symmetry in the bulk and on the supersymmetry breaking brane, 
the generated supersymmetry breaking masses can be flavor universal. The flavor symmetry is 
broken only on the other brane on which the Yukawa couplings are located. This setup leads 
to soft supersymmetry breaking masses that are completely general other than the fact that 
they are flavor universal. In particular, there is no imprint on the superparticle masses from the 
underlying gauge unification. 

It is, therefore, important to figure out the maximum value of one can obtain in generic 
high scale supersymmetry breaking scenarios. We first note that there is a "model-independent" 
source of the sensitivity of v to the fundamental parameters — the sensitivity of v to the gluino 
mass, M3, at M^nif. This is because M3 always gives contributions that make rrif grow at the 
infrared, which always pushes down the value of m|^^ at the weak scale. Since the sign of 
the contributions is definite, we cannot weaken this sensitivity by complicating renormalization 
group evolutions for a fixed value of M3. Now, there is a lower bound on M3 at Munif arising 
from the requirement that a sufficiently large At /mi is obtained at the weak scale, without 
making At at Munif extremely large and thus introducing a large sensitivity of v to At. We 
find M3(Miness) ~ 150 GcV, which leads to a factor of ~ 10 stronger fractional variation of v 
when we vary M3 at Munif: d\n.v/d\n.M^ k, 10. We conclude that, without a special relation(s) 
among various soft supersymmetry breaking parameters, the maximum value of in high 
scale supersymmetry breaking is 



^max 



10%. (7) 



This occurs in parameter regions in which \At/mi\ (9(1.5 ~ 2.5), ^ 250 GeV and Mg ^ 
450 GeV, where Mg is the gluino mass at the weak scale. In particular, the best points in Fig. 3 
both occur at At/mf ~ — 1.8 and = (^0, 

^2^ji/4 ^ 250 GeV. The electroweak VEVs satisfy 
tan /? ^ 5, and the Higgs boson mass is bounded by Mniggs ^ 120 GeV. 

If we want to improve fine-tuning further, we must consider a theory that gives smaller 
values of Mmcss, at least effectively, because an ultimate reason for the 10% tuning in Eq. (7) is 
the large logarithm In ( Mmcss /""^f) — In ( Munif /"^it). An important class of theories giving small 
Mmess is gauge mediation models [34, 35]. In these models, however, the size of A terms at 
Mniess is small, so that \At/mt\ ^ 1 at the weak scale. This requires large top squark masses to 
evade the Higgs boson mass bound, and thus leads to severe fine-tuning. In Fig. 4 we plot the 
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Figure 4: Contours of on the M^^ss-F / M-^^ess plane for the minimal gauge mediation models 
with TT-mess = 1 (left) and n^^ss = 4 (right) pairs of messenger fields in the 5 + 5* represen- 
tation. The instability bound implies the region in which the messenger fields are tachyonic, 
F/|MmessP > 1- The Higgs sector parameters are fixed as tan/3 = 15 and > 0. 



contours of A^^ as a function of F/Mmess and Mmess for /i > 0, in the minimal gauge mediation 
models with nmess = 1 and 4 pairs of messenger fields in the 5 + 5* representation of SU (5) D 
SU{?>)c X SU{2)l X U{1)y- Here, Mmess and F are the supersymmetric and supersymmetry 
breaking masses for the messenger fields, respectively. We find that fine-tuning cannot be better 
than 2% over the entire parameter region. 

2.3 Large At and small with low scale iVf mess 

Discussions in the previous subsections show that the fine-tuning problem in minimal supersym- 
metry is solved if soft supersymmetry breaking parameters are generated at a low scale with \At\ 
substantially (a factor of ~ 1.5 or so) larger than mj. It is, however, not so easy to achieve this 
in a simple manner. Suppose, for example, that the superparticles obtain masses through direct 
couplings to the sector that dynamically breaks supersjTiimetry at a scale A = 0(10 TeV). This 
will generate both A terms and non-holomorphic supersymmetry breaking squared masses, fh?, 
through operators of the form [[Z + Z^)$^$]5i4 and [Z^Z$''^$]04, where Z is the supersymmetry 
breaking (composite) chiral superfield and $ the quark, lepton and Higgs superfields. Flavor 
universality of the soft masses can be ensured by imposing a flavor symmetry on these couplings. 
Now, the strengths of these operators can be estimated using naive dimensional analysis [36]. We 
then find that the non-holomorphic masses are much larger than the A terms, im? (A/47r)^ and 
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A A/167r^, leading to an unwanted result of \At/mf\ l/An <^ 1. A possible way of avoiding 
this is to generate the operators given above in a perturbative regime at one and two loops, 
respectively, for example by using messenger-matter mixing in gauge mediation. This leads to 
Mmess in the 0(10 ~ 100 TeV) region, allowing fine-tuning to be relaxed to the (10~20)% level 
with a modest logarithm of ln(Mjncss/'"^t) — (3~6). 

From the viewpoint of obtaining a large A term at low energies, the simplest possibility is to 
have the operators Z"'')<l>'''<l>]g4 (and [Z'^ Z^'^^]0a) at tree level with 0(1) coefficients in units 

of the "cutoff" scale. Such a situation arises naturally if supersymmetry is broken associated 
with an extra dimcnsion(s) with the size of order (10 TeV)^^ [37 - 41] and if matter and/or 
Higgs fields propagate in the extra dimension(s). Consider that the MSSM gauge, quark and 
lepton superfields propagate in an extra dimension compactified on IZ2 with the length 7ri?. 
The two Higgs doublets are located on a brane, and the Yukawa couplings and the // term are 
introduced there. Then, if the boundary conditions for these fields are twisted by the SU{2)ji 
symmetry with an angle a, the theory just below 1/Ris the MSSM with the soft supersymmetry 
breaking parameters given by [40] 

a 2 f '^\^ A o'^ ^'c^ 

^1,2,3 = 1^Q,U,D,L,E = y-^) > ^«,<i,e = "2 — , (8) 

m'H^,H, = 0, I^B = 0, (9) 

where Mi 2,3 are the gaugino masses, and fnQ^u,D,L,E A^Ae the flavor universal squark and 
slepton masses and A terms, respectively. Taking a/ Rio be a few hundred GeV and 1/R — 
0(10 TeV), this gives a perfect spectrum for electroweak symmetry breaking: the messenger scale 
is low, Mjness — ^ / R, \At/mi\ is slightly smaller than 2 at the weak scale, and there is no strong 
hierarchy between the colored and non-colored superparticles. The origin of the small twist, 
a ~ (0.01 ~ 0.1), will lie in the dynamics of radius stabilization, as the SU {2)r twist in boundary 
conditions is equivalent to the supersymmetry breaking VEV in the radion supermultiplet [42]. 
A trade-off of this theory is that we lose a conventional picture of the supersymmetric desert, 
and thus a simple understanding of successful supersymmetric gauge coupling unification [43], 
although it might arise, for example, through some conformal property above 1/R with the 
conformality violation effect associated somehow with the zero-mode representations.^ In our 
view, the virtue of supersymmetric theories with a TeV extra dimension(s) lies in the fact that 
we can easily obtain large A terms at low scales. While the accommodation of the desert is 

^In such a scenario, the observed differences of the three low-energy gauge couplings should arise mainly from 
the differences of the gauge couplings in the bulk, and not from operators localized on a brane(s). Otherwise, the 
soft supersymmetry breaking parameters in Eq. (8) would receive large corrections from brane operators, and 
the colored superparticles would become much heavier than the non-colored ones, leading to the pattern which 
is not desirable in terms of electroweak symmetry breaking. 
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nontrivial, we think that these theories provide a much more natural solution to the so-called 
little hierarchy problem [44] compared with any other non-supersymmetric theories. 

It is possible to obtain a picture similar to the one described above without losing the con- 
ventional supersymmetric desert. In a theory where the moduli [45] and anomaly mediated [46] 
contributions to supersymmetry breaking are comparable [47 - 50], the mediation scale of super- 
symmetry breaking, Mmcss, can be effectively lowered without having a real physical threshold 
at Mjiicss [48]. The soft supersymmetry breaking parameters at M^^ss are then given essentially 
by those of boundary condition supersymmetry breaking or (equivalently) moduli mediated su- 
persymmetry breaking. (For a simple understanding of this property, see the Appendix.) One 
of the challenges to make a natural theory using this property is to have a sufficiently small fiB 
term to obtain a sufficiently large tan (3. Because of a large gravitino mass required to employ 
anomaly mediation, it is rather difficult to achieve the desired level of a (very) small /iB term. 

A model having all the desired features to have natural electroweak symmetry breaking 
keeping the supersymmetric desert has been given in Refs. [9, 10]. The effective messenger scale, 
-^mess, is lowered to the TeV region, a large At term with A-t/mi ~ —1.4 is obtained, and a 
large enough tan/3, tan/3 ^ 5, is accommodated by making small due to the renormalization 
group focusing effect and the elimination of the classical contribution. In the minimal setup, 
the soft supersymmetry breaking parameters similar to those of Eqs. (8, 9) are obtained at 
M^ess = 0(TeV): 

Mi,2,3 = Mo, ml^u,D,L,E = -f^ ^M,e = "Mq, (10) 

m]i^,H, = 0, liB^ 0, (11) 

where Mq is a parameter of order a few hundred GeV.^ Depending on the mechanism of // and 
^B term generation, the parameter B = ^iB/ ^ may have a non-zero value of order Mq/Att. As 
shown in Ref. [9], this theory does not suffer from fine-tuning, i.e. > 20%, and the region 
giving > 20% is consistent with the experimental constraints, such as the one from b — > 57, 
as long as the sign of fj, is positive [11].^ The essential point, again, is to generate a large \At/mf\ 
and moderately large tan /3 with a small (effective) messenger scale, Mmess ~ TeV. Note that 
small effective Mmess is achieved here by having special relations among various supersymmetry 
breaking parameters at Mimif, which allows us to evade the general result of Eq. (7). The top 
squarks should be light to eliminate fine-tuning, although they do not have to be as light as in 
the case of high scale supersymmetry breaking, because of rather small Mmess- A detailed study 
of the LHC signatures in this particular model will be given in section 4. 

^Note that the sign convention for the supersymmetry breaking masses adopted here, i.e. that of Ref. [29], is 
different from those in Refs. [9 11]. In particular, the sign of the A terms is opposite. 

^It is interesting to point out that the sign of is determined to be positive in the minimal setup, IJ.B = at 
Mmess, as long as Mmess is larger than the weak scale. 
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2.4 Characteristic spectra: light top squarks and light Higgsinos 

Considerations so far have highhghted certain generic features for a superparticle spectrum that 
leads to natural electroweak symmetry breaking in minimal supersymmetric theories. These 
have been obtained by considering mainly the tension between the m^Jj-g^ term in Eq. (3) and 
the Higgs boson mass bound from LEP II. Another important constraint on the spectrum comes 
from the term in Eq. (3). The fine-tuning arising from this term is about M^^^/2\ijl\'^, so 
that requiring A"^ > 20% (10%) leads to the bound ^ 190 GeV (270 GeV). This imphes 
that the Higgsinos should be much lighter than in typical mSUGRA or gauge mediation models, 
where relatively large top squark masses lead to a large /x parameter. 

Taking all these together, we find that a supersymmetric theory that has the minimal matter 
content and reproduces naturally the correct scale for electroweak symmetry breaking should 
have the following properties for the superparticle spectrum: 

• The A term for the top squarks is large, \At/mf\ (1.5 ~ 2.5) at the weak scale. This leads 
to a large mass splitting between the two top squarks: 

—rut^ (1.5 ~ 2.5)777,4. (12) 

Here, we have assumed m^^ ^ ''T^u-i- "^^^ splitting as small as ~ 772,4, however, may be 
allowed if Mmcss is small, Mmcss = O(TeV). 

• The top squarks should be light, i.e. = {rn,Q,^mfj_J^^^ should be small, to reduce the 
sensitivity of v toyt. How small mi should be depends on the value of Mmcss and the amount 
of A^^ required. For the case of high scale supersymmetry breaking with A~^ pa 10%, the 
bound is very strong, m^- ^ 300 GeV, leading to rri^^ ~ 100 GeV. 

• Small values for the top squark masses imply that we cannot push up the Higgs boson mass 
much larger than the tree-level value. Typically, we find 

Mffiggs ^ 120 GeV. (13) 

• The ratio of the electroweak VEVs should also be moderately large 

tan/? ^5, (14) 

to have a sufficiently large Higgs boson mass at tree level. This imphes that the /jB term 
should be (significantly) smaller than 2|//p -|- 171^^ + mjj^ at the weak scale. 

• The fjL parameter should be small 

< 190 GeV (270 GeV), (15) 

for A~^ > 20% (10%). This leads to light Higgsinos, which may be the LSP, or if not, may 
significantly mix with the LSP. 
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Figure 5: Characteristic spectra for the superparticles which give the correct scale for electroweak 
symmetry breaking without significant fine-tuning. The spectrum in (a) arises typically in a high 
scale supersymmetry breaking scenario with ^ 10%, while that in (b) arises in a theory 
where supersymmetry is broken by boundary conditions, or moduli contributions, with small 
(effective) Mmess- 



The features described above still leave several possible patterns for the superparticle spec- 
trum, which can lead to somewhat different situations. In Fig. 5 we depict possible patterns 
which are representative for generic cases. In Fig. 5(a), we depict a pattern that typically arises 
in a high scale supersymmetry breaking scenario with A^^ ^ 10%, for example in the mSUGRA 
model with m|^^,m|^^ ^ ml (see Fig. 3). In this situation the top squarks are rather light, 
TUf = ("^Q^m^g)^/^ ^ 300 GeV, with the light top squark close to its experimental bound, 
m^^ ^ 100 GeV. The Higgsinos, h, are also light, although ^ 10% allows as large as 
~ 270 GeV. In principle there are little constraints on the other gaugino and sfermion masses, 
except that the gluino should be heavier than about 450 GeV and that the squarks cannot be 
much lighter than the gluino. The universal gaugino mass relation, M^/ = M2/ g"^ = Mi/gf, 
may or may not be satisfied. The LSP will be either bino-like, wino-like, Higgsino-like, or a 
mixture of these states. If the LSP consists mainly of the bino, its thermal relic abundance can 
give the correct dark matter density either through a large mixture with the neutral Higgsinos 
and/or wino, coannihilation with the stau or stop, or resonant annihilation through s-channel 
Higgs boson exchange. Note that some of these options are not available in typical mSUGRA 
points, in which top squarks are heavy and thus the fi parameter is large. For other cases of 
wino-like and Higgsino-like LSPs, the production should be nonthermal. It is interesting that a 
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large Af term also makes it easier to satisfy the constraints from the precision electroweak data 
with small top squark masses [51]. It is important to perform detailed LHC studies of this class 
of spectra, although a large ambiguity for the gaugino and sfermion masses will make a thorough 
study of the parameter space somewhat complicated. These spectra will also be interesting for 
a future e+e~ linear collider [52]. 

If we want to reduce fine-tuning further, for example to eliminate it altogether (A^^ > 20%), 
we must generate the soft supersymmetry breaking parameters at low energies. Because of a 
small logarithm, ln(Mjncss/'"^?)) constraints on the top squark sector is slightly weaker in this case: 
lAt/m^l ^ 0(1) and rrif ^ (9(700 GeV). On the other hand, the constraint on the /i parameter is 
stronger because of a stronger requirement on A^^: for A^^ > 20%, we obtain the bound on the 
Higgsino masses rrif^ ^ 190 GeV. In fact, generating these spectra with small Mmess is nontrivial, 
and one of the ways is to adopt (effectively) the scheme of boundary condition, or moduli, 
supersymmetry breaking at low energies. This generically leads to a rather ordered spectrum 
at low energy, e.g. universal gaugino and sfermion masses at the weak scale. This situation is 
depicted in Fig. 5(b), where almost degenerate gaugino and sfermion masses are assumed. The 
relation between the gaugino and sfermion masses is model dependent. An important implication 
of this class of spectra is that the LSP is one of the neutral Higgsinos, unless the gravitino is 
fighter. To identify this LSP to be the dark matter, it must be produced nonthermally, for 
example, as in [53]. As discussed in [11], such dark matter can have an interesting implication 
on direct dark matter detection experiments sucfi as CDMS II. 

In tfie rest of tfie paper, we focus on studying LHC signatures for the latter class of spectra, 
given in Fig. 5(b). We do so partly because it gives milder (or no) fine-tuning, and partly because 
there is little LHC study directly related to this case. In particular, we first focus on signals 
expected from the Higgsino LSP at the LHC, and discuss under what conditions the signals are 
most useful. We then demonstrate in section 4 that the signals can indeed be used in realistic 
analyses, using the explicit model discussed at the end of subsection 2.3. The determination of 
model parameters are also discussed there, which may be useful to discriminate between various 
possible models. 

3 Higgsino LSP at the LHC 

We have seen that in a large class of theories where the supersymmetric fine-tuning is solved, 
the Higgsinos are the fightest among the superpartners of the standard model particles. A 
characteristic pattern for the superparticle spectrum in these theories is depicted in Fig. 5(b). 
The three gauginos are almost degenerate at the weak scale, as well as the squarks and sleptons. 
The ratio of the gaugino and the sfermion masses is model dependent. In this section, we identify 
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characteristic signals of these spectra at the LHC, and discuss under what conditions the signals 
are useful in realistic analyses. We assume throughout that the gravitino is not lighter than the 
Higgsinos, so that the LSP is the lightest neutral Higgsino.^ 



3.1 Dilepton invariant mass distribution from X2 



Xi I'^l in the 



Higgsino LSP scenario 

An important feature of the spectra depicted in Fig. 5(b) is that there are three almost degenerate 
neutrahno/chargino states, Xi: X2 Xi: with the masses ^ ^ 190 GeV. If the gaugino 
masses are sufficiently larger than |//|, which we assume to be the case, these states are almost 
purely the Higgsinos, with the mass splittings given by 



where we have assumed a moderately large tan l3, e.g. tan l3 ^ 5, 6w is the Weinberg angle, and 
Ml and M2 are the U{1)y and SU{2)l gaugino mass parameters, respectively. This implies that 
X2 and Xi undergo three-body decays to x5- particular, X2 has the leptonic decay mode 



At hadron colliders, this decay mode can give important information on the properties of the 
initial- and final-state ncutralinos [3, 54 - 56]. Below we show that the dilepton arising from 
the decay of Eq. (18) can provide an important test for the Higgsino nature of the lightest two 
neutralinos and X2- 

The three-body decay X2 ~^ Xi occurs through the diagrams shown in Fig. 6. In the limit 
where the mass difference between the two neutralinos is much smaller than the Z boson mass 
and where the slepton masses are much larger than the decaying neutralino, m^o — m^o <^ niz 
and mi^^mi^ 3> ?7t.^o, the effects of these diagrams are described by a single low-energy XiX2^"'"^~ 
four-Fermi operator for each chirahty of leptons. This implies that the distribution shape of the 
dilepton invariant mass, Mu = a/ {pi+ -\-Pi-Y, is completely determined by the masses of the 
two neutralinos, Xi and X2) where pi+ and pi- are the four- momenta of Z+ and . Let us now 
ado})t tlie i^liase eom eiitioii in liic li all the mass eigenvalues and the mixing matrix elements 

^In fact, none of our analyses changes unless the gravitino mass, TO3/2, is smaller than 0(1~ 10 keV) because 
the lightest Higgsino then lives long enough so that it can be treated as a stable particle for collider purposes. 
For a smaller gravitino mass, the lightest Higgsino decays into a Higgs boson and a gravitino, followed by the 
Higgs boson decay h — > hb. This can be used to measure the Higgs boson mass, e.g., by selecting four 6-jet events 
and plotting M^^ invariant masses. For 7713/2 = 0(0.01 ~1 keV), we may also have displaced bb vertex signals. 




(16) 



(17) 




(18) 
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Figure 6: The diagrams contributing to the ~^ Xi^~^^ decay. 



for the neutrahnos are taken to be real. This basis can always be taken as long as there is no CP 
violating effect in the neutralino mass matrix, which we assume throughout. We then obtain the 
following dilepton invariant mass distribution after performing appropriate phase space integrals: 



dMn 



oc M, 



^ K§ - "^iO + K? + ^W^fi - 2Mti + Qr^^m^om^oMti ,(19) 



for < Mil < ~ ^^'^ dT{y^ x^l^l^) / dMu = for Mu > m^o — m^o. Here, m^o = 
I M^o I and m^o = \ M^o \ are the absolute values for the two smallest neutralino mass eigenvalues 
M^o and M^o with m^o < m^o, and r/^ = sgn(M-j^o) sgn(M^o) is the relative sign between them. 
In fact, with the LEP II bound on the Higgsino masses, the assumption of m^o — m^o <^ mz 
implies that the two neutralinos are nearly degenerate: Am = m^o — m^o <^ ^xl- The Mu 
distribution is then further simplifies to 

"^^^^"rflrf ^^^'^ ^ {2^^' + + 3^x)M/}, (20) 

for 0<Mii< Am and dT{xl Xi l^l')/dMii = for Mu > Am. 

There are two important features for the Mu distribution in Eq. (20) which can be used to 
test the Higgsino LSP scenario. First, the endpoint of the distribution, m^o — m^o, is expected 
to be very small: 

2 

M-- = m^o - m^o ^ ^ = 0(10 GeV), (21) 

where we have set Mq = Mi ~ M2. Given the LEP II bound on the chargino mass, such a small 
mass splitting between Xi ^"^^ X2 cannot arise in a theory where the LSP is gaugino-like and 
the three gauginos respect the universal mass relation, M^/gl = M2/g1 = Mi/gf. Second, the 
Higgsino LSP necessarily leads to the opposite signs between M^o and M^o, so that = — 1 in 
Eqs. (19, 20). This is because in the gauge eigenbasis the 2x2 neutral Higgsino mass matrix 
takes a purely off-diagonal form ((0, — /i), (— /i, 0)), which gives one positive and one negative 
eigenvalues after diagonalization including the effects of mixing with the gaugino states. The 
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Figure 7: Dilepton invariant mass distribution from the X2 ~^ Xi decay. Solid lines represent 
the distributions for the Higgsino LSP case (left panel) and the case of gaugino-like x5 X2 
with sgn(Mi) = sgn(M2) (right panel). The curves with dashed and dotted lines represent the 
ones obtained from Eq. (20) with 77^ = —1 and +1, respectively (in both panels). 

resulting distribution is quite different from the one with 77^ = +1, which arises in the case of 
the gaugino-hke LSP with sgn(Mi) = sgn(M2). 

In Fig. 7 we plot Mu distributions calculated for several different choices for the soft super- 
symmetry breaking parameters. In the left panel, we plot the Mu distribution for the Higgsino 
LSP case by a solid curve. The curve is drawn using the complete expression with tan/3 = 10, 
Ml = M2 = 480 GeV and /i = 170 GeV, which gives m^o — m^o ^ 20 GeV. We have varied 
the slepton masses mj = mj^ = mj^ in the range 1.001 < rj = mj/m^o < 10, but it does not 
lead to any visible change of the curve. In the figure wc have also drawn curves obtained using 
the approximate expression of Eq. (20) for rj^ = —1 (dashed line) and 1]^ = +1 (dotted line). 
We find that the expression of Eq. (20) with rj^ = —1 well approximates the full result. The 
small discrepancy arises from corrections higher order in [Am/mz)'^ ■ The right panel shows the 
Mil distribution in the case of gaugino-hke Xi ^-^d X2 with sgn(Mi) = sgn(M2) (solid lines). 
The parameters are chosen to be tan/3 — 10 and /i — 500 GeV, and Mi and M2 are chosen 
such that the same values of m^o and m^o as in the left panel are obtained: Mi — 158 GeV 
and M2 — 182 GeV. The slepton masses are varied as rf = rrif/m^o — 10, 1.1, 1.01, 1.001, and 
we find that the Mu distribution in this case depends on but only when it is very close to 
the X2 mass. As in the left panel, we also draw curves obtained from Eq. (20) with rf^ — —1 
(dashed line) and r]-^ — (dotted line). We find that the approximate curve with r]^ — well 
reproduces the full result with r^- ^ 1.1. 
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The plots in Fig. 7 clearly show that the Higgsino LSP case can be discriminated from the case 
with Tj^ = sgn(M^o) sgn(M^o) = +1, even with the same mass difference m-^o — m^o, regardless 
of the values for the other supersymmetry breaking parameters such as the slepton masses. In 
particular, we find that the behavior of the Mu distribution near the endpoint is completely 
different between the two cases of r^^ = —1 and +1. This can be understood in terms of the 
selection rule for the orbital angular momentum due to the CP properties of the two neutralinos 
Xi and X2 [57]- For r]^ = +1 (—1), the Mu distribution near the kincmatical endpoint, which 
corresponds to the limit of slow moving Xi in the X2 ^^st frame, should give an S-wme (P-wave) 
behavior, leading to Mu oc {Am — MuY^"^ {Mu oc {Am — MuY^"^) near the endpoint. While 
gaugino-like and X2 could potentially give a similar distribution if sgn(Mi) = — sgn(M2), the 
shape of the Mu distribution together with the smallness of the endpoint can provide a powerful 
tool to test the Higgsino LSP scenario considered here, which necessarily leads to Am <S mz 
and r]^ — —1.^ 

To demonstrate that the signatures discussed above are really useful at the LHC, we must 
check that there are no other leptons from the superparticle cascade decays which bury the 
signatures. We must also show that the shape of the Mu distribution from the x^ — > x° 
decay is preserved under selection cuts in the analysis at a level that different shapes for 77^ = +1 
and —1 can be discriminated in a realistic detector. We address the first issue in the next 
subsection, in the context of a class of theories discussed in subsections 2.3 and 2.4. The second 
issue will be addressed in section 4, where we explicitly demonstrate that the Mu distribution 
can indeed be used to test the Higgsino LSP scenario, as well as to extract the information on 
the neutralino masses, by performing Monte Carlo simulations using a specific theory. 

3.2 Higgsino LSP with quasi-degenerate gauginos and sfermions 

Let us consider the pattern of the superparticle masses depicted in Fig. 5(b), which can naturally 
arise in a theory where moduli-type, or boundary condition, supersymmetry breaking is employed 
with small Mmess- Specifically, we consider the following spectrum. The three gauginos are almost 
degenerate at the weak scale with the masses denoted by mg. The squarks and sleptons are also 
nearly degenerate with the masses denoted by m^, although the two top squarks, ti and 
can have substantially different masses because of a large mass splitting due to large Ai. The 
mass splittings among the gauginos and among different squarks and sleptons arise at higher 

^It is interesting to point out that the signatures discussed here arise only from the fact that the neutral 
Higgsinos are pseudo-Dirac fermions. The same technique, therefore, can also be used to test the idea of pseudo- 
Dirac gauginos, depending on the spectrum for the other superparticles (for examples of theories giving pseudo- 
Dirac gauginos, see [58]). For instance, if the LSP is one of the pseudo-Dirac bino or wino states, similar 
signatures may arise in the dilepton invariant mass distribution, depending on the existence of other sources of 
leptons and/or patterns of cascade decays for the superparticles. 
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order, but they are expected to be small and of 0(10%). The /i parameter is smaller than about 
200 GeV, so that it is smaller than both nig and rriq. 

We first consider the case with rriq < rrig. In this case, a squark cannot decay into a gluino, 
so that it decays as g — > XiQ.' , XiQ x^q. Here, q and q' represent quarks having the same and 
different flavors with g, and the lightest chargino, Xi i ^"^^ fhe lightest two neutralinos, x52' 
the charged and neutral Higgsinos, respectively, with small mixings with the gaugino states. On 
the other hand, the gluino, once produced, decays into a quark and a squark g qq, followed by 
squark decay discussed above. In these decay chains, leptons arise only from decays of xt and X2'- 
Xi — ^ Xi aiid X2 ~^ Xi ■ In the Mu distribution analysis, we select events having two and 
only two leptons with the same flavor and opposite charge. The number of dileptons arising from 

decays is then small over the relevant energy region of Mu = 0(10 GeV), compared with that 
from X2 decay. (The background from xt decays can actually be estimated using opposite-sign 
opposite-flavor leptons and thus subtracted using the combination e^e" + jJi^ — ^~ — /x+e".) 
The only remaining issue then is the squark branching ratio into x^ and the x^ branching ratio 
into leptons. In the parameter region we consider, the former is typically of 0(10%) and the 
latter is ~ 3% for both e^e~~ and modes, which are large enough to produce an appreciable 

number of dilepton events.^ We therefore conclude that the signatures discussed in the previous 
section can be used at the LHC for rriq < rrig. This claim will be confirmed in the next section, 
where we perform an explicit study using Monte Carlo simulations. 

In the case of rriq > rrig, a squark decays mainly into a quark and a gluino, although it may 
also decay into a wino or bino by a small amount. The gluino then undergoes three-body decays: 
g — >■ xtqq' 1 XiQQ or X2qq. The branching ratios for g Wqq^'\ Bqq may also be non-zero if 
the gluino is slightly heavier than the wino and/or bino, but these modes are highly suppressed 
by smallness of the phase space. The gluino branching ratio into X2 is of 0(10%), and the x^ 
branching ratio into leptons (including both e and yu) is ~ 6%, implying that an appriciable 
number of dilepton events can be obtained from xH decay. The background from Xi decay is, 
again, not important. The only sources of leptons that could potentially destroy the signatures 
are decays of W and B, produced by squark decay. The dominant decay modes of W and B, 
however, will be into a Higgsino (either one of xt 1 Xi and xV) ^nd a W, Z or Higgs boson, so 
that the dangerous modes giving I'^l" directly are suppressed by the three-body phase space and 
a small gaugino-Higgsino mixing. The number of dangerous dileptons from W and B decays, 
therefore, is at most of the same order as the ones from x^ decay. Since the Mu endpoint for these 

^Strictly speaking, opposite-sign same-flavor dileptons may also arise from the three-body decay of the gluino 
g Wqq followed by the wino decay W ^ II ^ HXi^ if tlic gluino is slightly heavier than the wino, e.g. by about 
0(10%), due to higher order effects. The branching ratio of this mode, however, is extremely small because of a 
large phase space suppression, so that the resulting dileptons are completely negligible compared with the ones 
arising from X2 decay. A similar comment also applies to q ^ qll followed by I — > Ixi- 
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dileptons is about rrig — \iJ^\ and the distribution is suppressed for Mu much smaller than this 
values, it is unlikely that these leptons destroy the signatures from X2 decay. (We select events 
having two and only two leptons when doing the Mu analysis.) We thus find that the Higgsino 
LSP signatures discussed in the previous subsection are also useful in the case of rriq > rrig at 
the LHC. 

We finally consider the case where rriq = rrig at the leading order (at tree level). In this case 
we expect that the masses of the gluino and squarks, rrig and rrig, are slightly (0(10%)) larger 
than those of the wino, bino and sleptons, m^, and mj, due to higher order (radiative) 
effects. The orderings among rrig and rriq and among m^, and mj are model-dependent. 
With these spectra, gluinos and squarks once produced decay mostly into W or B plus a few 
jets, although a small fraction decays directly into light Higgsino states, xt^ Xi and X2- Decays 
of the electroweak gauginos differ depending on the ordering of m^r, rrij^ and m^. If there is a 
slepton with the mass smaller than that of a gaugino, e.g. < m^f, there is a sizable branching 
ratio for the gaugino decaying into the slepton, W ^ II. This gives a large amount of opposite- 
sign same-flavor dileptons through the slepton decay I — > Zx°, which can potentially destroy the 
signatures from x^ decay. On the other hand, if all the sleptons are heavier than W and B, 
these gauginos decay mainly into a Higgsino (one of x^, and X2) ^'^'^ ^ ^ or Higgs boson. 
The branching ratio into X2 is typically of 0(10%), and the desired signatures are obtained from 
X2 decay. There are other dileptons from three-body decays of W and B, such as W ^ Xi^h 
but the number of these dileptons is sufficiently small. We thus expect that the Higgsino LSP 
signatures are useful for rriq — rrig as long as all the sleptons are heavier than the electroweak 
gauginos. 

We have seen that the Higgsino LSP signatures discussed in subsection 3.1 are useful, i.e. 
not buried by dileptons from other superparticle decays, in a large class of theories motivated by 
solving the supersymmetric fine-tuning problem. In the next section, we explicitly demonstrate 
that these dilepton signatures can indeed be used in realistic analyses by performing Monte 
Carlo simulations in a theory with rrig > rriq. We also present a technique which can essentially 
determine all the superparticle masses in a class of theories discussed here, up to small mass 
splittings of 0(10%) among different squarks and sleptons and a smaller splitting between the 
electroweak gaugino masses. 

4 Natural Super symmetry at the LHC 

In this section, we perform a Monte Carlo study for a class of theories discussed in the previous 
section and in subsections 2.3 and 2.4, which naturally leads to the correct scale for electroweak 
symmetry breaking. We explicitly demonstrate that the signatures discussed in the previous 
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section can be used to test the Higgsino LSP and extract the small mass difference between the 
two neutral Higgsinos at the LHC. We also devise a series of cuts that allows us to determine 
all the relevant superparticle masses in theories with rrig > rriq: m^o, rriq, Am and rrig. In 
particular, we perform the analysis in the context of the model based on mixed moduli and 
anomaly mediated supersymmetry breaking, discussed at the end of subsection 2.3, and show 
that the model can be tested at the LHC, up to theoretical uncertainties of a; 15% on various 
superparticle masses. 



4.1 Framework 



The basic setup for our analysis is the same as that in subsection 3.2 (and subsection 2.3). The 
three gauginos are almost degenerate at the weak scale, rrig = Mi ~ M2 — M-^. and the squarks 



and sleptons are also nearly degenerate, rriq = {mnY^'^ ~ (m^)-*^/^ ^ (m|,)-'^/^ ~ (ml)-*^/^ 



(m 



E 



. The A parameters are nearly universal at the weak scale, A = A„ ~ ~ Ag, with 
A satisfying |A/mg| ^ 0(1). The (top) squark masses, mg, should not be very large, and the 
ratio of the electroweak VEVs should satisfy tan /3 ^ 5. In the analysis in this section, we only 
consider the case rrig > rriq. 

To perform an explicit Monte Carlo study, we must choose particular parameter points. For 
this purpose, we take the model based on mixed moduli and anomaly mediated supersymmetry 
breaking, discussed at the end of subsection 2.3, and choose the parameters within the region 
satisfying the condition > 20%. In particular, we take [9, 10] 

( A'P \ 

Mi,2,3 = Mo, ml^^^^^^E = A^Ae = -Mo, m\^^^^ - O (^^J , (22) 

at some scale Mmess = 0(100 GeV~TeV), where 

450 GeV ^ Mo ^ 900 GeV, (23) 

for tan /3 ^ 20. For smaller tan /5, the lower bound in Eq. (23) increases; for tan /3 = 10 (5) the 
lower bound becomes ?a 550 GeV (900 GeV). In our analysis, we take /x, txia and tan/3 to be 
free parameters in the Higgs sector, which are left undetermined after the electroweak symmetry 
breaking condition is imposed on /x, /^S, m|^^ and rr{\^^. As shown in Ref. [11] the constraint 
from the b ^ sj decay chooses the sign of /i to be positive. We thus have 

sgn(/x) = +l, |/x|<190GeV, ^ 300 GeV, tan/3^5. (24) 

The last three conditions come from the naturalness criterion, > 20%. 

In our Monte Carlo study, we choose two parameter points given in Table 1, satisfying 
Eqs. (22, 23, 24). The point I is representative for the case with a relatively low superparticle 
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point I 


point II 


Mo [GeV] 


600 


900 


H [GeV] 


170 


170 


niA [GcV] 


250 


250 


tan (3 


15 


15 


-^mess 


1 TeV 


1 TeV 



Table 1: Two representative parameter points of the model used in Monte Carlo simulations. 

mass scale, while the point II for the case with a high superparticle mass scale. These points 
satisfy the experimental constraints such as the ones coming from the Higgs boson mass and the 
b ^ s'-f decay, within theoretical uncertainties. We note, however, that the constraints from the 
Higgs boson mass and 6 — > 57 are not very important in our present context, because they are 
sensitive to the parameters in the Higgs sector, such as rriA and tan (3, whose precise values are 
not relevant in our LHC study below. 

The physical masses for the superparticles are obtained from the inputs in Table 1 as follows. 
We interpret the input masses as the running masses in the DR' scheme at the scale Mmess, and 
evolve them down to the superparticle mass scale using renormalization group equations. We 
then add the effects of electroweak symmetry breaking, such as the D-term contributions to the 
scalar masses and the gaugino-Higgsino mixings, and convert the running masses to the pole 
masses by including finite threshold corrections, using the code SuSpect 2.3 [28]. The resulting 
superparticle masses are given in Table 2 for the two parameter points in Table 1. Strictly 
speaking, this procedure is not quite meaningful in the context of the model under study, because 
we generically expect unknown higher order corrections of 0(l/87r^) in the expression of Eq. (22), 
which can be comparable to some of the low energy corrections included here. Nevertheless, this 
procedure allows us to incorporate the fact that the colored particles are systematically heavier 
than the non-colored ones by 0(10%), which we expect to hold in realistic situations. We thus 
perform our Monte Carlo study using the masses given in Table 2, although one should remember 
that there are intrinsic theoretical uncertainties of 0(10%) for the superparticle masses in the 
model. 

To perform the analysis, we generate both supersymmetric and standard model events us- 
ing Pythia 6.324 [59]. We generate supersymmetric events for the two parameter points in 
Table 1. The number of events generated for each point is equivalent to the integrated lumi- 
nosity of 30 fb""^, which corresponds to the three-year running of the LHC at low luminosity. 
The superparticle decays are calculated using the code SDECAY 1 . la [60], with the results trans- 
ferred to Pythia and used in the event generation. For the estimation of standard model 
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Table 2: Superparticle masses in GeV for points I and II in Table 1. The masses for the second 
generation squarks and sleptons are not listed because they are nearly degenerate with the 
corresponding first generation squarks and sleptons. 

background, we have generated 0.5M QCD 2— >2 events for each bin of the transverse momen- 
tum: 100 GeV <Vt < 200 GeV, 200 GeV < Vt < 400 GeV, 400 GeV < Vt < 800 GeV, and 
Pt > 800 GeV. We have also generated the l^+jets events with W — > ei/, //i/, rv (0.5M events 
for 50 GeV < Vt < 200 GeV and 0.2M events for pr > 200 GeV), the Z+jets events with 
Z iyp,T+T- (0.5M events for 50 GeV < pr < 150 GeV and 0.2M events for pr > 150 GeV), 
IM events for the tt production, and 0.2M events for each ZZ, ZW and WW production. 
These standard model events are simply scaled to 30 fb"^^ when estimating standard model 
backgrounds. While our background estimations arc correct probably only to a factor of a few 
due to inherent uncertainties associated with the QCD effects, we expect that our analysis is 
not much affected by this because the standard model background can be pretty much reduced 
by our cut selections, as we will see later. Some of the analysis, e.g. the effective mass analysis 
in subsections 4.4 and 4.5, may be affected by these uncertainties, but then we can always raise 
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the cut on E^^^^ and recover the usefulness of the analysis. 

For the detector simulation, we use AcerDET 1 .0 [61], a generic fast detector simulation and 
reconstruction package for the LHC, which has a similar principle of operation to the official fast 
simulation package of the ATLAS detector, ATLFAST [62]. The package performs identification 
and isolation of leptons, photons and jets in terms of detector coordinates: pseudorapidity 1], 
azimuthal angle 0, and cone size AR = \J + (Zir;)^. Lepton, photon and jet four-momenta 
arc smeared, and the cluster selections arc made based onpj- and \!]\. Isolation criteria arc applied 
to leptons and photons in terms of the distance from other clusters, AR > 0.4, and of maximum 
transverse energy deposited in cells in a cone AR = 0.2 around the cluster. The calibration of 
jet four-momenta is also performed, and each jet is labeled either as a light jet, 6-jet, c-jet or 
r-jet, using information from event generators. (We use default parameters for these selection, 
isolation, calibration and labeling processes.) For the 6-jet identification, we further implement 
6-tagging efficiency of 60% per a 6-labeled jet, with mistagging probability of 10% for a c-labeled 
jet and 1% for a fight jet. For the r-jets, we use efficiency of 50% per a r-labeled jet, with the 
mistagging probabifity of 10% for other jets. 

For each event, we apply the following trigger selections [62]: one isolated electron with 
Pt > 20 GeV, one isolated photon with > 40 GeV, two isolated electrons/photons with 
Pt > 15 GeV, one muon with pt > 20 GeV, two muons with pt > 6 GeV, one isolated electron 
with Pt > 15 GeV and one isolated muon with px > 6 GeV, one jet with px > 180 GeV, three 
jets with Pt > 75 GeV, and four jets with pt > 55 GeV, where isolated electrons/photons, 
isolated muons and jets must be in the central regions of pseudorapidity \r]\ < 2.5, 2.4, and 3.2, 
respectively. In our analysis, we consider only events passing one of these criteria. 

In our study, we ignore possible systematic errors caused by cuts and choices of fitting 
functions and regions, and take only into account the statistical errors. Based on good agreements 
between the input values and the fit results obtained in subsections 4.3, 4.4 and 4.5, however, 
we expect that these neglected errors are not much larger than the statistical errors included in 
the analysis. 

4.2 Characteristic features: short cascades and fewer leptons 

We begin by identifying characteristic features relevant to the LHC arising for the superparticle 
spectra under consideration. At the LHC, most of the superparticle productions comes for 
rnqjUig ^ 1 TeV from squark and gluino productions. For point I (point II) in Table 1, the 
cross sections are ~ 3.8 pb, 25.3 pb and 23.0 pb (0.27 pb, 2.98 pb and 3.71 pb) for the gg, 
gq and qq productions, respectively, where the total superparticle production cross section is 
about 55.2 pb (9.22 pb). We find that the squark-gluino pair production and the squark pair 
production dominate in our parameter space. After the production, a squark decays mostly 



26 



X? 

Figure 8: Decay cascades used in the analysis. 



xS 



into xf or Xi and a quark, q — > Xil'^ Xi?; but it also decays into ^ xtly with a small 
branching ratio of 0(10%). The xt ^.nd X2 produced then decay into Xi and quarks/leptons 
through three-body decays, giving leptons with branching fractions of 0(10%). For the gluino, 
it decays into a quark and a squark with the branching ratio of ~ 100%, which is followed by 
squark decay. 

An important feature of these decays is that the decay chains are relatively short. Compared 
with the case where the wino decay into a slepton is open, for example, decay chains with the 
present spectra are shorter in average. Another important feature is that the number of leptons 
arising in the cascades is significantly smaller than in the case where colored and non-colored 
superparticles have a large mass hierarchy, e.g. as in a typical mSUGRA parameter point. 
We can thus use (# of 1 lepton events)/(# of lepton events) to make a first guess that the 
superparticles may have a spectrum like the one considered here. In fact, the features described 
here allow us to test certain generic aspects of the spectra, such as the nature of the LSP, and to 
determine the basic mass parameters in a simple manner. These information can then be used 
to test or discriminate between possible models, as will be discussed in subsection 4.6. 

In our analysis, the following decay cascades are used: q — > X2q x^l~^l~q, q XiQy 
and g qq Xilli which are depicted in Fig. 8. Here, we have not discriminated between 
quarks and antiquarks. Using the kinematics of these cascades, we can determine the masses 
of the gluino, squarks and neutralinos, as well as the small mass difference between x\ aiid X2-> 
model independently. This will be shown in subsections 4.3 and 4.4 for the case of point I and 
in subsection 4.5 for the case of point II. Various kinematical endpoints, such as the ones for 
dileptons, two leptons plus a jet, and a combination of two jets, will be used. Precisions of order 
a few to ten percent are achieved, as will be shown later below. 
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4.3 Determination of the squark and neutral Higgsino masses 

In this subsection we show that using the kinematics of the cascade decay q — > ~^ Xi^~^^~Q 
and those for the squark pair production with q — > XiQj we can determine rriq, m^o and m^o 
without using input from particular models. This information can thus be used for nontriv- 
ial tests of the model predictions, as will be discussed in subsection 4.6. The analysis also 
demonstrates that the Mu distribution discussed in section 3 is indeed useful in testing the 
Higgsino nature of the lightest neutralinos. In this and the next subsections, we use point I 
in Table 1 (Mq = 600 GeV) for the analysis. The same analysis will be repeated for point II 
(Mo = 900 GeV) in subsection 4.5. 

We first look at the Mu distribution from the three-body decay X2 ~^ Xi^^^ - discussed 
in section 3 the endpoint and the shape of the distribution measure the mass difference of the 
two neutralinos, Am, as well as the relative CP property of the two neutralinos, 77^^. 

We select the dilepton events with the following cuts: 

• Ef^^ > 300 GeV 

• At least two jets with pt > 50 GeV 

• Two and only two leptons with the same flavor and opposite charge 

• Veto 6-jets 

In Fig. 9, we show the Mu distribution obtained with the cuts described above. The standard 
model backgrounds are effectively reduced by the cuts (hatched histogram). We can clearly see 
the endpoint of Mu around 15 GeV, which can be the first test for the Higgsino LSP scenario. 
Note that while Mu is small, Pt^s of the leptons are not so small because of the parent X2^s 
transverse momenta, so that these leptons are not much affected by the trigger selections of 
subsection 4.1. 

We have performed a fit of the distribution by using the theoretical curve including the effect 
of the finite Z-boson mass, with a linear background distribution. In this fitting process, the 
standard model background events are smeared in order not to artificially magnify the statistical 
uncertainty due to the scahng of the generated events to 30 fb~^. A reasonable fit is possible 
only for the r]^ = —1 case (solid line), and we obtain the endpoint value 

M;^^ = 15.30 ± 0.15 GeV, (25) 

which is consistent with the input value of the mass difference Am = 15.40 GeV (see Eq. (21)). 
The theoretical curve with the rj^ = +1 case is superimposed in the plot (dashed line). We 
can clearly see the deviation from the simulated distribution, especially near the endpoint. We 
thus conclude that the smallness of the endpoint, M^^^, together with the shape of the Mu 
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Figure 9: The distribution of the dilepton invariant mass, Mu, for point I. Hatched histogram 
represents the standard model background, which is smeared over five bins in order not to 
magnify the statistical uncertainty due to the scaling of the events. The solid line is the best fit 
function for the signal plus background, obtained using the theoretical curve with rj^ = —1. It 
can be clearly distinguished from the r/^ = +1 case, drawn by the dashed line. The endpoint is 
extracted to be 15.30 ± 0.15 GeV. 

distribution characteristic of rj^ = —1, provide an extremely powerful test for the Higgsino LSP 
scenario. If these signals are actually observed, they strongly suggest that the LSP is one of the 
nearly degenerate neutral Higgsinos. 

At this stage, we only have information on the neutralino mass difference, but further infor- 
mation can be obtained by stepping up the cascade, i.e. by combining dileptons with the quark 
jet from the squark decay. We can construct two independent Lorentz-invariant quantities Mug 
and Miq, whose endpoint values are given by 

[^-^J ■ (26) 

The two endpoints coincide because the final state leptons come from the three-body decay. The 
measurement of the two endpoints, therefore, can give only one additional information. 

For the event selection for the Mug and Mig measurements, we have imposed a cut on Mu 

• Mu < 15 GeV 

in addition to the cuts for the Mu measurement, in order to reject incorrect lepton pairs. The 
jet which is combined with the lepton(s) is selected from the two largest pt jets. We choose the 
one that gives the smaller Mug to see the endpoint of the distribution. 
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Figure 10: The distributions of Mug (left) and Mig (right) for point I. Hatched histogram repre- 
sents the standard model background, which for Mug is smeared over five bins in order not to 
magnify the statistical uncertainty due to the scaling of the events. The solid line is the best fit 
function for the signal plus background near the endpoint, obtained by using a linear function 
with Gaussian smearing for the signal and a linear function for the background. The endpoint 
is extracted to be 179.6 ±5.8 GeV for Mug. The Mig endpoint cannot be extracted clearly using 
a simple linear function. 

In Fig. 10, we show the distributions of Mug (left) and Mig (right). The endpoint is clearer in 
Miiq than in Mig. We thus use the Mug endpoint for the mass determination. We have performed 
a fit of the Mug distribution near the endpoint, with a linear function with Gaussian smearing. 
A linear function is also assumed for the background distribution. The best fit function is drawn 
in the figure, and we obtain 

M,7/" = 179.6 ± 5.8 GeV. (27) 

This is consistent with the expected endpoint obtained using the formula of Eq. (26), which is 
~ 180 GeV. We have also tried to estimate the endpoint in the Mig distribution. A reasonable fit, 
however, cannot be obtained with linear functions. Better functions are needed if one wants to 
extract the endpoint from the Mig distribution. From the Mn and Mug endpoint measurements, 
we now have two relations among three mass parameters, rriy^o, m^o and rrig. We still need one 
more independent quantity to determine all the three masses. 

In principle, the threshold value of Mug with the cut Mu > ^M^^^ (0 < < 1) could provide 



30 



the required additional information: 
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However, with the hmited statistics and the narrow physical Mu region of the Higgsino LSP 
scenario, the threshold of Mu^ is not quite useful for the mass determination. With fixed 
and M^^^, ^Ifq^ has a very little sensitivity to the mass parameters. 

We therefore have to look for another quantity to determine the three masses. Such a quantity 
can be obtained by analyzing the squark pair production process followed by the two squarks 
decaying into two jets and two x?'s. Although we cannot reconstruct the squark four- momenta 
due to two escaping invisible neutralinos by the event by event analysis, we can extract a relation 
between rriq and m^o by the endpoint analysis of the Mt2 variable defined in Ref. [63]. This 
variable is designed to take the maximal value at the squark mass when we input the correct 
m^o in the calculation. The definition is given by 

M|2= . min . [max{m|(p^^p-^),m|(p^^p??^-)}], (29) 

--.miss I „miss --.miss 

where p^ and p^^ are the transverse momenta of the jets from the squark decays, and p™'^^ is 
the missing transverse momentum. The transverse mass, m|,, is defined by 

^r(PT, Pr) ^ml + ml + 2{E^E!^ - p« • p^). (30) 

By identifying the endpoints of Mt2 for various input values of m^o, we can obtain a relation 
between rrig and m^o, which can provide the last information to determine the three masses, 
m^o, m^o and rriq. 

To select the squark pair production events, we use the following cuts: 

• Ef^^ > 300 GeV 

• Veto leptons, 6-jets, r-jets 

• Two and only two jets with px > 50 GeV 

For our assumptions on the 6-tagging and r-tagging efficiencies, see subsection 4.1. 

In Fig. 11, we show the Mj'2 distribution for the input value of m^o = 200 GeV as an example. 
We can see a clear edge in the distribution around 500 GeV. Fitting with a linear function with 
a linear background, wc obtain the endpoint 505.3 ± 0.6 GeV. With the rich statistics (63,859 
events survive the cuts) and the sharp edge, we can measure the endpoint quite accurately. Note 



31 




Figure 11: The Mj'2 (defined in the text) distribution for the input value of the neutrahno 
mass m^o = 200 GeV for point I. Hatched histogram is the standard model background. The 
endpoint is extracted by fitting the signal plus background histogram with a linear function, 
and the background near the endpoint by a linear function. The endpoint is obtained to be 
505.3 ± 0.6 GeV. 

that this is not the squark mass itself — it is the value of some quantity that would become the 
squark mass if our hypothetical input neutralino mass is in fact the true neutralino mass. 

Combining the MJ^^^ and M™^^ measurements with the Mt2 analysis, we can now determine 
all the three masses, m^, m^o and ^ti^o, as follows. First, the endpoint analysis of the cascade 
decay, M™^^ and M^^^, gives two constraints on the three mass parameters, leaving one param- 
eter unfixed. If we take this parameter to be rn^o^, we can draw a curve on the m^o-rrig plane, 
using the constraints from M™*^^ and M^^*^^. On the other hand, the Mj'2 analysis of the squark 
pair production gives another relation between m^o and m^, giving an independent curve on the 
same plane. The intersection of the two curves will then give the real values of [m^o, rriq). In 
Fig. 12 we show the two curves explained above with the statistical errors. The Mt2 curve is 
obtained by performing the Mt2 endpoint measurements for six different input values of m^o 
and then interpolating them with a smooth curve. The measured values of m^o and rriq by this 
combined analysis are indicated by shaded bands with the la statistical errors. 

The effectiveness and accuracy of the method are demonstrated in Fig. 13. To draw the 
figure, we have generated 10,000 "experiments" and considered that in these experiments the 
values of M™^^, Mj]^^^ and M^^^ are determined according to the Gaussian distributions with 
the statistical errors given in Eqs. (25, 27) and by the Mt2 fit. We have then calculated the 
three mass parameters using the method described above and have plotted their distributions. 
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Figure 12: Two curves on the m^o-rriq plane deduced from the cascade decay analysis, M™^^ 
and M]^^, and the squark pair production analysis, M^2^, for point I. Both curves are obtained 
by inputting hypothetical values of fn^o, which is taken as the horizontal axis. The intersection 
determines the real values of m^o and rrig. The obtained masses with the la statistical errors 
are shown by shaded bands. 

These plots show that m^o and rriq have larger tails in large mass regions. This represents the 
fact that the two curves in Fig. 12 are more similar in a larger m^o region than in a lower region. 
The input values of the mass parameters are indicated in Fig. 13 by arrows, and we find that the 
correct values are obtained within reasonable statistical uncertainties. By fitting the histograms 
with the Gaussian distribution, we obtain 

m^o = 169 ± 17 GeV, m^- = 486 ± 11 GeV, Z\m = 15.30 ± 0.15 GeV. (31) 

This demonstrates that the neutralino and squark masses can be measured with 10% and 2% 
level accuracy, respectively, at the LHC in the Higgsino LSP scenario. The mass difference 
between the two neutral Higgsinos can be measured at 1% accuracy. The information on these 
masses are very useful to test particular models, as will be discussed in subsection 4.6. 

4.4 Determination of the gluino mass 

With the knowledge of the squark and neutralino masses, we can determine the gluino mass 
using the kinematics of the g ^ qq x\qq cascade decay. The invariant mass of the two final 
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Figure 13: The statistical errors of the measured mass parameters m^o, rUq and Am for point I. 
Arrows indicate the input values used in the Monte Carlo event generation. The accuracy at 
the level of 10%, 2% and 1% is obtained for m^o, rUq and Am, respectively. 
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Therefore, if the endpoint of the Mjj distribution arising from this cascade is measured, we can 
determine the gluino mass. 

In supersymmetric models, the gluino production is mostly from the ^ + g or ^ + ^ pair 
production, which necessarily gives additional jets from the other side of the squark/gluino 
decay. Those additional jets cause an uncertainty for the selection of the correct jet pair. In 
order to reduce this combinatorial background, we select the events with three hard jets, which is 
typical in the ^f + g production, and choose a pair of jets which gives the smallest Mjj among three 
combinations such that the calculated Mjj would not exceed the endpoint in Eq. (32). (Note 
that this docs not mean that the selected pair is necessarily the correct one, but it guarantees 
that the endpoint of the Mjj distribution is given by the right formula, Eq. (32).) 

The selection cuts we use are the following: 

• Ef^^ > 300 GeV 

• Veto leptons, 6-jets, r-jets 

• Three and only three jets with > 50 GeV 

With these cuts, we obtain the Mjj distribution shown in the left panel of Fig. 14. The endpoint 
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Figure 14: The Mjj (left) and M^h (right) invariant mass distributions for point I. In the Mjj 
analysis, we select events with three hard jets and take a combination that gives the smallest 
Mjj, out of the three combinations, in order to see the endpoint. A fit is performed with a linear 
function with Gaussian smearing together with a linear function for the background events. The 
endpoint is obtained as 377 ± 10 GeV. In the M^b analysis, we select events with three hard jets 
including two hard 6-jets. A similar endpoint, 370 ± 11 GeV, is obtained using the same fitting 
function. 

structure is visible around 400 GeV. By fitting the histogram near the endpoint by a linear 
function with Gaussian smearing and a linear background shape, we obtain 

Mj^^" = 377 ± 10 GeV. (33) 

The large value of shown in the plot is caused by the artificially magnified statistical uncer- 
tainty of the standard model background due to the scaling of the events. We have checked that 
the reasonable value of is obtained without the standard model background. 

We can perform the same analysis for the g bb ^ Xi^b decay by requiring two 6-jets. We 
show the resulting M^h distribution in the right panel of Fig. 14. In this case, we do not suffer 
from the combinatorial background, although the statistics is reduced. If we assume ~ rriq, 
which is indeed the case here, we can use this endpoint for the gluino mass reconstruction. We, 
however, do not use this analysis in the following because it will make our analysis more model 
dependent. For the purpose of testing the particular model, however, the Mbb analysis can be 
used as a consistency check (or to extract some information on the value of tan/3). 

Combining the information of Mjj'^^ with the analysis in the last subsection, we can determine 
the gluino mass. The reconstructed gluino mass is shown in Fig. 15. The input value of rrig = 
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Figure 15: The statistical error of the gluino mass for point I. The arrow indicates the input 
value used in the Monte Carlo event generation. The accuracy at the level of 2% is obtained. 

623 GeV is indicated by the arrow. We see that the reconstruction is successful within the 
statistical uncertainty. With the approximation of the statistical fluctuation to take the Gaussian 
form, we obtain 

m-g = 632 ± 13 GeV. (34) 

We find that a quite accurate (fa 2%) measurement of the gluino mass is possible by this method. 

There is another method of extracting the gluino mass, which can be used in any model 
within the class considered here. This is to use the effective mass Mes defined by 

Meff = £;r' + ^p^, (35) 

i 

where the sum runs over all the jets. The peak location of this variable is known to have a 
correlation with the gluino and squark masses [2, 3]. In particular, as we will sec below, wc have 
a definite relation between the peak location of Meg and the superparticle masses within the 
model used here. To perform this analysis, wc use the cut criteria listed in Ref. [64] to select the 
events, except for the lepton, 6-jet and r-jet vetoes: 

• > 4 jets with pr > 50 GeV 

• > 2 jets with pr > 100 GeV 

• Ef^^ > maxjlOO GeV, 0.25 

• Transverse sphericity St > 0.2 
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Figure 16: An empirical relation between the peak location of Mgfj and a combination rrig + 
rrig — m^o. We have generated 50,000 events for each 30 sample points in the mixed moduli- 
anomaly mediation model. A very good linear relation is obtained. Note that this relation will 
be modified if different selection cuts are used. The standard model background has not been 
included in searching the peak location. 

. Z^0(p^,p2) < 170° 

• Veto leptons, 6-jets, r-jets 

By generating supersymmetric events for various parameter points in the model, we find an 
excellent linear relation between M^g^^ and nig + rrig — m^o as shown in Fig. 16. The relation is 
given by 

M^g^^ = 285 GeV + 0.763 {rrig + m-g - m^o). (36) 

Note that this should be regarded as a sort of theoretical prediction as we have not included the 
standard model background in producing the plot. By using this empirical fact, we can extract 
the combination rrig + rrig — m^o from the M^s peak analysis. 

In Fig. 17 we show the distribution of Mes described above. We obtain the peak location 
M^^^^ = 977 GeV by fitting the histogram near the peak with a Gaussian function. With the 
assumption that the theoretical relation in Eq. (36) is accurate at a 5% level, we obtain the 
gluino mass 

m-g = 590 ± 62 GeV, (37) 

by combining the M^^^^ analysis here with the analysis in the last subsection. We find that the 
error amounts to 0(10%). We thus conclude that the Mjj endpoint analysis is more useful than 
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Figure 17: The distribution of Meg for point I. The QCD background is smeared in order not 
to artificially magnify the statistical uncertainty due to the scaling of the events. The location 
of the peak is determined by fitting with a Gaussian function near the peak. It is given by 

the M^^^"^ analysis to determine the gluino mass not only because it is more model independent 
but also because it has better accuracy. 

4.5 The case with large superparticle masses 

In this subsection we repeat the analyses in subsections 4.3 and 4.4 for the case of large super- 
particle masses (point II in Table 1) to see if accurate measurements are still possible despite 
the smaller statistics due to smaller superparticle production cross sections. The analyses below 
show that essentially the same method can be used to determine the mass parameters in good 
accuracy. In fact, similar or even better accuracy is obtained for the rrtq and rrig determination 
compared to the case with low superparticle masses, as we will see below. 

The Mil and Muq distributions are shown in Fig. 18. We have used the same cuts with those 
in subsection 4.3 for the event selection. We can see the clear endpoints in both distributions. 
The Mil distribution is fitted with the theoretical curve with 77^ = —1. A good fit is obtained 
only for rj^ = —1, allowing a successful measurement of 7]^. The endpoints are obtained as: 

^max ^ g 4g ^ g 21 Qgv, Ml^^"" = 223 ± 12 GeV. (38) 

where the Muq endpoint is determined by a fit using a linear function with Gaussian smearing 
together with a linear background shape. 
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Figure 18: The distributions of Mu (left) and Mug (right) for point II. Hatched histogram 
represents the standard model background, which is smeared over five bins in order not to 
magnify the statistical uncertainty due to the scaling of the events. For Mu, we have searched 
the endpoint by fitting with the theoretical curve assuming rj^ = —1. The endpoint is obtained 
to be 9.48 ± 0.21 GeV. For Mug, we have used a linear function with Gaussian smearing for the 
signal and a linear function for the background. The endpoint is obtained to be 223 ± 12 GeV. 

The Mt2 distribution with an input m^o = 200 GeV is shown in Fig. 19. The signal to 
background ratio is not so large, but there is no significant background near the endpoint because 
of the large squark mass. The endpoint measurement, therefore, does not suffer seriously from 
the standard model background. We obtain 716. 6± 1.9 GeV for this value of the input neutralino 
mass. We then repeat the analysis for six different input neutralino masses and obtain a curve 
on the m^o-rriq plane by interpolating these points. 

The two curves obtained from the Mu and Mug endpoints and the endpoint are shown 
in Fig. 20. We find that the crossing angle of the curves is larger compared to the case in Fig. 12. 
This property makes it possible to measure the squark mass in this point with similar accuracy 
to the case of low superparticle masses, even though we have larger statistical uncertainty. 

For the Mjj endpoint measurement, we use a different cut for the jet pt- We take 

• Three and only three jets with pt > 100 GeV 

instead of px > 50 GeV, used in the analysis of the low superparticle mass case, because a 
clear endpoint is not obtained with the 50 GeV cut. The resulting Mjj distribution is shown in 
Fig. 21. By fitting with a linear function with Gaussian smearing and a linear background, we 
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Figure 19: The distribution for the input value of the neutrahno mass m^o = 200 GeV for 



point II. Hatched histogram is the standard model background. The endpoint is extracted by 
fitting the signal plus background histogram with a linear function, and the background near 
the endpoint by a linear function. The endpoint is obtained to be 716.6 ±1.9 GeV. 
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Figure 20: Two curves on the m^o-rriq plane deduced from the cascade decay analysis, 
M™^^, and the squark pair production analysis, M^2^, for point II. Both curves are obtained 
by inputting hypothetical values of m^o . The intersection determines the real values of m^o and 
iTLq. The obtained masses with the Icr statistical errors are shown by shaded bands. 
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Figure 21: The Mjj invariant mass distribution for point II. We select events with three hard 
jets and take a combination that gives the smallest Mjj, out of the three combinations, in order 
to see the endpoint. A fit is performed with a linear function with Gaussian smearing together 
with a linear function for the background events. The endpoint is obtained as 609 ±13 GeV. 

obtain 

Mj^^" = 609 ± 13 GeV. (39) 

Combining all the results, we can determine the four mass parameters, m-^o, rriq, Am and rrig. 
The estimation of the statistical uncertainties is given in Fig. 22. The input values are indicated 
with the arrows, which are all within reasonable statistical uncertainties. The gluino mass is 
obtained with slightly larger values. This is caused by the systematics that the Mjj endpoint 
tends to give larger values than the one obtained in Eq. (32) when we use the pr > 100 GeV 
cut for jets. Therefore, the estimation of the systematic error will be important in this analysis. 
By fitting the histograms with Gaussian functions, we obtain 

m^o = 164 ± 24 GeV, m^- = 700 ± 12 GeV, Z\m = 9.5 ± 0.2 GeV, (40) 

m-g = 940 ± 15 GeV. (41) 

The neutralino, squark and gluino masses can be measured at 15%, 2% and 2% level accu- 
racy, respectively. Since this point represents the case of the highest superparticle masses from 
the naturalness requirement (see Eq. (23)), the above analysis shows that the method of mass 
determination developed here covers the entire region of the parameter space. 

For completeness, we show in Fig. 23 the M^s distribution. The standard model background 
is huge, although it is not so significant around the peak location. By using the relation in 
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Figure 22: The statistical errors of the measured mass parameters rn^o, rriq, Am and rrig for 
point II. Arrows indicate the input values used in the Monte Carlo event generation. The 
accuracy at the level of 15%, 2%, 2% and 2% is obtained for m-^o, m^, Am and m^, respectively. 

Eq. (36), we obtain the gluino mass 

= 801 ± 85 GeV, (42) 

which is slightly deviated from the input value of 917 GcV (by lAa). This little discrepancy is 
mainly caused by the shift of Mge in the lower direction due to the standard model background. 
For a realistic use of the Mcs analysis in the high superparticle mass region, one needs to develop 
a better understanding of the shape of the standard model background and/or to devise a better 
cut (especially on E^^^^) to reduce the standard model background. Note that the relation in 
Eq. (36) will be modified if different cuts are used. 

4.6 Testing the model with mixed moduh-anomaly mediation 

One of the most important features of the LHC experiment is its potential of ruling out models. 
With the limited precisions of various measurements, it is extremely important to develop meth- 
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Figure 23: The M^s distribution for point II. The QCD background is smeared in order not 
to artificially magnify the statistical uncertainty due to the scaling of the events. The peak 
location is determined by fitting with a Gaussian function near the peak. It is found to be 
MPg^'^ = 1306 GeV. 

ods of testing model predictions rather than just measuring parameters under the assumption 
of a particular model. We here take the model with mixed moduli-anomaly mediation as an 
example and discuss possible ways of testing the model. 

We have already seen that one of the characteristic features, the Higgsino LSP, can be tested 
using the distribution of the dilepton invariant mass from X2 decay. Other interesting features of 
the model include approximate universality of the gaugino masses and the definite ratio between 
the sfermion and gaugino masses, given in Eq. (22). We can test these features by using the 
mass parameters obtained in the previous analyses. Specifically, we can calculate the parameter 
Mq in three different ways and compare them with each other. If the relations predicted in the 
model hold, the three values must coincide within (mostly theoretical) uncertainties/corrections. 

The first way of calculating Mq is to use the measured neutralino mass difference. Am. By 
using Eq. (16) with Mi ~ M2 ~ Mq, we can extract Mq as 

(43) 



Am 

This must give the correct value of Mq up to corrections of ~ 15%, which come mainly from 
0{\n\ sin2/9/Mo) corrections in diagonalizing the neutralino mass matrix and from the effect of 
running between the effective messenger scale, Mmess, and the gaugino mass scale. The other 
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Figure 24: A nontrivial test for the mixed moduli-anomaly mediation model in the case of 
Mo = 600 GeV (left) and Mq = 900 GeV (right). Three ways of calculating Mq, i.e. ml /Am, 
'\/2mg and m^, give the same value within ^ 15% theoretical uncertainties. For m^, we have 
plotted the values obtained in two different ways, i.e. from MJ^"'' and M^^''^. 

two ways use direct measurements of the squark and gluino masses: 

Mo ~ V2mg, Mo ~ mg, (44) 

which must also give the correct value of Mo up to corrections. The corrections come, for 
example, from running between Mmess and the superparticle mass scale as well as from finite 
supersymmetric QCD corrections. For the squark masses, the SU{2)l and f/(l)y -D-terms also 
give corrections of O(m|/Mo). All these corrections, again, amount to ~ 15%. We thus conclude 
that the three values of Mo calculated using Eqs. (43, 44) must coincide within 15%, if the model 
is actually realized. This can provide a rather nontrivial test for the model. 

We show in Fig. 24 the values of Mo calculated in three different ways for point I (Mo = 
600 GeV; left panel) and for point II (Mo = 900 GeV; right panel). The shaded regions indicate 
the ±15% range around the true values of Mo (= 600 GeV and 900 GeV). For mg, we have 
plotted both values obtained from the Mjj endpoint and the Mes peak analyses. As we can see, 
all measurements agree with each other within 15% undertainties, both for the cases of point I 
and point II. Moreover, we can even understand a nature of the dispersions in the plot. We find 
that the 0{\fi\ sin2/5/Mo) correction to m|/Z\m is negative whereas the QCD corrections to the 
squark and gluino masses are positive, as expected from theory. Since experimental errors on 
determining these quantities are rather small, the structure of these dispersions might be useful 
for deducing further detailed structures of the underlying theory, such as the structure of higher 
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order corrections to the superparticle masses. 

5 Discussion and Conclusions 

Weak scale supersymmetry is often said to be fine-tuned, especially if the matter content is 
minimal. Is it true? If the LEP II bound on the Higgs boson mass pushed the top squark masses 
above a TeV, as is sometimes stated in literature, this would be a true statement. The size of 
the top-stop loop contribution to the Higgs boson mass-squared parameter, ni'^H^, would then 
be larger than about (250 GeV)^ even for a unit logarithm, ln(Mnicss/'"^t) ~ 1, leading to fine- 
tuning of order 10% independently o/an underlying mechanism of supersymmetry breaking (see 
Eqs. (2, 3)). In fact, in most supersymmetry breaking models, the logarithm is (much) larger 
than 1, leading to (much) severer fine-tuning. For example, fine-tuning is already as bad as 2% 
for Mmess ~ 100 TeV. Such heavy top squarks, therefore, would not allow natural electroweak 
symmetry breaking in the context of minimal supersymmetry. 

The situation, however, is not as described above if there is a large At term at the weak 
scale. For \At/mt\ ^ (1.5~2.5), the top squark masses can be as small as mf ^ (250~400) GeV 
to evade the LEP II constraint of Mniggs ^ 114.4 GeV, so that the "model- independent" con- 
tribution to rnjj^ from top-stop loop is only about (100 GeV)^, even including the contribution 
from the At term. (The lower bound on here, ^ 250 GeV, arises in fact from the direct 
search bound of the top squark, ^ 100 GeV, and not from the Higgs boson mass bound.) 
Such light top squarks allow "fine-tuning" to be only of 0(20%) for Mmess = 0(10 TeV) and 
allow completely natural electroweak symmetry breaking for smaller values of Mmess- With these 
low energy spectra, the amount of fine-tuning, A^^, can also be small for theories with large 
logarithms. We find that A^^ ^ 10% is possible even in theories with Mmess ~ Munif- There- 
fore, under the current experimental constraints, the supersymmetric fine-tuning problem should 
not be regarded as the problem of minimal supersymmetry itself but as the problem of specific 
supersymmetry breaking mechanisms. 

Characteristic features required to have (relatively) natural electroweak symmetry breaking in 
minimal supersymmetric theories are (i) a large A term for the top squarks, \At/mt\ ~ (1.5 ~ 2.5) 
(ii) light top squarks (iii) a moderately large ratio of the electroweak VEVs, tan /? ^ 5, and (iv) 
a small n term, < 190 GeV (270 GeV) for A-^ > 20% (10%). A generic implication of these 
low energy spectra is a relatively fight Higgs boson, Mniggs ^ 120 GeV. There will be classes 
of theories leading to these features/spectra, and we have identified two representative ones 
(see Fig. 5). Among them, a class giving nearly universal gaugino and sfermion masses at low 
energies (see Fig. 5(b)) can make electroweak symmetry breaking most natural. Examples for 
these theories are obtained by employing moduli-type, or boundary condition, supersymmetry 
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breaking (effectively) at a low scale. A consistency with the desert can be explicitly recovered if 
the setup is realized through mixed moduli and anomaly mediated supersymmetry breaking. 

An important consequence of the class of theories described above is that the Higgsinos 
are the lightest among the superpartners of the standard model particles. Assuming that the 
gravitino is not much lighter than the Higgsinos, which is actually the case in the model with 
mixed moduli-anomaly mediation, the existence of the nearly degenerate Higgsino states (x?, 
X2 and xt) '^^^ interesting signals at the LHC The signals arise in the invariant mass 
distribution of dileptons arising from X2 decay: a smallness of the endpoint and a particular 
shape determined by the relative CP property of the two neutralinos, Xi a-^d X2- We have 
argued that these signals are indeed useful in a wide variety of circumstances within the class of 
theories considered here. 

We have explicitly demonstrated the usefulness of the signals in realistic analyses by per- 
forming Monte Carlo simulations, including detector simulations and background estimations. 
We have also presented a method that allows the determination of all the relevant superpar- 
ticle masses, m^o, niq, rrig and Am — m^o — m^o, independently of the details of the model. 
This allows us to determine all the superparticle masses within the class of models considered, 
up to theoretical uncertainties of ~ 15%. Note that some of the existing techniques, e.g. the 
Miiq threshold analysis, is not very useful here because of the near degeneracy of the Higgsino 
states. We can, nevertheless, determine the four mass parameters with the precisions of order 
a few to ten percent, by combining various endpoint analyses. This is extremely important 
because it provides ways to test various possible models, which generically give nontrivial rela- 
tions among these parameters. We have demonstrated this in the case of the model with mixed 
moduli-anomaly mediation, and shown that the model can indeed be tested (and thus can be 
discriminated from models that give significantly different relations among the four parameters). 

It will be possible to perform further tests for the class of models discussed here. An important 
issue is to measure the A parameters, especially that of the top squarks. This may be done, for 
example, along the lines presented in [65]. There is also an important interplay between collider 
physics and cosmology. As discussed in Ref. [11], the present class of models has a large discovery 
potential in ongoing and future direct dark matter detection experiments, such as CDMS H, if 
the Higgsino LSP composes the dark matter of the universe. Now suppose that the mass and 
the detection cross section for the dark matter are measured in (one of) these experiments. The 
results from the LHC can then be used to perform a consistency check on the LSP mass and to 
provide a constraint on the Higgs sector parameters, and tan /3, because the detection cross 
section depends strongly on these parameters. Together with the other data from the LHC, 
such as the one for the Higgs boson mass, we will be able to determine all the parameters of the 
model with certain accuracy. 
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In performing all these analyses, inputs from a particular model(s) will be very impor- 
tant/useful, especially if one wants to pin down the parameter point of minimal weak scale 
supersymmetry. The model(s) assumed is then better to be a "likely" one, compared with other 
possible models that can also "accommodate" the same set of measurements. Naturalness of 
electroweak symmetry breaking, together with the simplicity of a model, will then keep playing 
an important guiding role in these "model selection" processes, which will most likely take the 
form of processes of "slowly convincing ourselves." 
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A Renormalization Group Properties of Moduli Media- 
tion Models 



We here study interesting renormalization group (RG) properties of moduli mediated supersym- 
metry breaking models. It is found in Ref. [48] that mixed moduli-anomaly mediation models 
have an interesting RG property — the contributions from anomaly mediation can be canceled 
by the actual one-loop running effect at a certain energy scale, mimicking a pure moduli me- 
diated model at that energy scale. Since a clear derivation of the result is not available in the 
paper, we present it here. We show that the result can be understood as a rather straightforward 
consequence of special properties of moduli mediated models. We use the method of "analytic 
continuation into superspace" [66], which provides a powerful way of analyzing RG eqiiations 
in softly broken supersymmetric theories. We explain the basic mechanism of this interesting 
property and the origin of the conditions for the cancellation to happen. We also point out that 
the effects of mass thresholds are under a good theoretical control. 

We consider the case where the effective Lagrangian is given at a scale A as follows: 




(45) 
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where T is a spurion field with a non- vanishing F-component. This is the form of the La- 
grangian obtained in moduh (radion) mediated models. The lowest component of T represents 
the volume of the extra dimensions in which the gauge fields propagate. The rational number rj 
represents how much fraction of the extra dimensions the matter Qi propagates, compared with 
the gauge fields. For example, if the gauge fields and the matter field Qi propagate in six and 
five dimensional spacctimc, respectively, rj is given by (5 — 4)/(6 — 4) = 1/2. 

We now exploit the following property of moduli mediated models to show certain special 
RG properties of these models. We first note that at tree level there are following simple scaling 
relations associated with the rescaling of the moduli field T: 

S^aS, Zi^a^'^Zi for T ^ aT, (46) 

where S and Zi are the gauge kinetic function and the wavefunction factor defined hj C 3 
[SW^Wfilr and C 3 [ZiQ\Q-^£,, respectively. Wc then find that these scaling relations can be 
extended to the one-loop level if the moduli rescaling, T aT, is supplemented by the following 
rescalings of the RG scale /ir and the Yukawa couplings: 

In ^ ^ a In ^, A,,, ^ a^'-^+^^+^-i)/^ A,,„ (47) 

where Xijk are the superpotential Yukawa couplings appearing in Eq. (45) (the "physical" Yukawa 
couplings are given by yijk = Xijk/ {ZiZjZkY/"^). Once this property is proved (see later), we 
can use these scaling relations to show that the gauge kinetic function, S, and the wavefunction 
factor, Zj, take the following form: 



(48) 



Z. = (T + Ttr..Z.(^^^JM_,_2_,„^), (49) 

at an arbitrary scale /xr. Here, S and Zi are some functions, and we have used the fact that 
Zi can depend only on T through the combination T + because of the invariance of the 
Lagrangian under the transformation T ^ T -\- 

Now, suppose that the condition + + = 1 is satisfied for the fields having the Yukawa 
interaction Ajj^. We then find that the T dependencies in S and Zi appear only with the 
renormalization scale ^r. In this case, the gaugino masses, A terms, and soft scalar squared 
masses are simply given by: 



mx = ^[In-Sjf = Mo 



2hg' ,^ ^R 



In — 

(47r)2 A 



(50) 
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1 + 2(7i + 7,- + 7fe) In 



A 



m2 = -[lnZi]B = Mo2 



r. + 47.1n^ + 27.(ln/^)^ 



(51) 



(52) 



at an arbitrary renormalization scale ^r. Here, g and fe represent the gauge couplings and the 
beta function coefficients, d\n{l/ g'^)/d\n.^ji = — 26/(47r)^, respectively, and 7j and 7j are the 
anomalous dimensions, dlnZi/dlnfiR = — 27j, and their derivatives with respect to the scale 
I^R, ii = dji/dln iiji. The overall supersymmetry breaking parameter Mq is defined by 



Mo 



(53) 



where the subscript A denotes the lowest component. Note that the soft supersymmetry breaking 
parameters in Eqs. (50 - 52) are given by simple functions of the quantities at the scale /ir, the 
gauge couplings, beta functions and anomalous dimensions, as well as ln(/Xij/A). 

It is rather simple to prove the scaling relations in Eqs. (46, 47). The RG equations for the 
gauge couphngs (the gauge kinetic functions S — l/2g'^ -\ ) are given by 



where we have defined t = Ir{ij,r/A). This obviously leads to 

^S{T,t) = -f—aS{T,t) = --^ 
dt ^ ' d{at) ^ ' (47r)2 

On the other hand, the RG equations for SiaT^ at) are given by 



(54) 



(55) 



d 



d{at) 



S{aT, at) 



(47r)^ 



(56) 



We thus find that aS{T, t) satisfies the same RG equation as S{aT, at). With the initial condition 
S{T, 0) = T, we can determine the integration constant: 



S{aT,at) = aS{T,t). 



(57) 



This imphes that S scales as S ^ aS for {T,t) — > {aT,at), which is the relation we wanted to 
prove. In fact, the result here is not a special property of moduli mediated models. The effective 
Lagrangian in softly broken supersymmetric theories can always be recast in the form of Eq. (45) 
as far as the gauge sector is concerned, so that the expression for the gaugino masses in Eq. (50) 
is a (well-known) general result. 
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We can prove the scaling properties of Zj's along the same lines. The RG equations for Zj's 

are 

I In = -27, = \ Y: - 2(^ + Sr'C^\R\ (58) 

j,k * ^ ^ 

at one loop, where C2\R) is the quadratic Casimir operator for the superfield Qj.^ We can 
again transform this to 

d , , ^ ^ ^ 1 a^^l A,- 



j,k •' 

= -Yt—77T7 r. [/ r.. -2(s(aT,at) + S(aT,aty) ci'\R),(59) 



where we have used Eq. (57) in the second equation. This equation shows that a'^^Zi{T, Xijk,t) 
satisfies the same RG equation as Zi{aT,a^^^~^^^~^^''~^^/'^Xijk,at). With the initial condition of 
Zi{T, Xijk, 0) = [T -\-T'^Y\ the integration constant is determined and we obtain 

Zi{aT, a(^^+^^+^'=-i)/^A,,fe, at) = oJ-^Zi{T, A,,fc, t). (60) 

This is the scaling relation for Z^ given in Eqs. (46, 47). 

There is a subtlety if the gauge group contains a C/(l) factor. In this case the Fayet-Iliopoulos 
term 

'(i^e^VY, (61) 



may be induced at one loop, which contributes to the soft scalar squared masses as 

ml^-[hiZ,]D + gUY,, (62) 

where gy is the U{1) gauge coupling and is the U{1) charge of the superfield Qi. The RG 
equation for ^ is given by [67] 



d -2 
Je ~ (4^ 



E^^^i-^^^ + H^E^- (63) 



Since the combination of '^■YilnZi is RG invariant, i.e. i^7i = 0, ^ is never generated if 
• Yi In Zi = Q (and ^ = 0) at the classical level. Therefore, we can neglect the contributions to 
mf from the Fayet-Ihopoulos term if the condition ^ . Fjrj = is satisfied. 

^Precisely speaking, S + in Eq. (58) should be the real gauge coupling superfield R defined by -R — 
(Tg/Stt^) \nR = S + S^ - Ei(2"i/87r2) InZj (in the NSVZ scheme). The difference, however, is irrelevant at the 
one-loop level. 
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In the case where there is a mass threshold M, we find that the scahng properties of 
Eqs. (46, 47) are maintained if the rescahngs of T, /ir and Xijk are supplemented by 

, M , M 

In— ^aln— . 64 
A A 

The functions S and Z^, in this case, can depend on ln(M/A)/T and ln(|M|7|A|2)/(r + Tt), 
respectively. The soft supersymmetry breaking terms, therefore, obtain additional contributions: 

Am.^^Moln-, (65) 

M 

AAijk = -2(Z\7i + + ^7fe)Mo In — , (66) 



Am', = Ml 



4A^iln — + 2A^i Mn — 



(67) 



where Ab, A'ji and Aji are the changes of b, ji and 7^ at the scale M ((high scale value) — 
(low scale value)). The gauge coupling g in Eq. (65) is the one at the scale fiR. 

The derivation here should make it clear the origins of the special properties of Eqs. (50 - 52) 
and the required condition + rj + = 1. Since the RG equations for the gauge and Yukawa 
couplings take the form of dg/dt ~ g^ and dy/dt ~ + yg"^ at one loop, and g'^ oc 1/T and 
y"^ oc {ZiZjZk)^^ oc (T + T^)"*^^*"'"'"^"'"''*^ in moduli mediated models, it is clear that one-loop RG 
equations are invariant under the rescaling (T, t) — > (aT, at) if + + = 1 is chosen. This 
scaling property then guarantees the forms of Eqs. (48, 49), leading to Eqs. (50 - 52). (This 
also makes it clear that these properties persist under the existence of arbitrary generational 
mixings.) This simple scaling property clearly cannot persist at higher loop orders, so that the 
properties of Eqs. (50 - 52) are that of one-loop RG equations. 

Inclusion of anomaly mediation is straightforward at this point. We should simply replace A 
in Eqs. (48, 49) by A$, where $ (= 1 + 771.3/26'^) is the chiral compensator field. (Note that the 
compensator field $ does not couple to T as T is a dimensionless chiral superfield.) A curious 
similarity between moduli and anomaly mediations is manifest here. In particular, at the scale 

InA^!!!^, (68) 
m 2Mo' ^ ' 

the F-component of Indiu/ A^) /T as well as F- and £)-components of ln(//|/|A$p)/(T' + T'^) 
vanish. Therefore, either S ox Z does not have F- or £)-components if the condition ri-\-rj-\-rk — 1 
is satisfied. The solutions of the RG equations at this scale are remarkably simple: 

mx = Mo, Ajk = -Mo, ml=r,Ml (69) 

If there is a mass threshold, the solutions are obtained by simply adding the contributions in 
Eqs. (65 - 67), because of the ultraviolet insensitivity of anomaly mediated contributions. 
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